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ABSTRACT 


(ii) 


The effect of diffusion limitation on selectivity of 
some complex reactions in a cylindrical catalyst pore has beer 
studied using Langmuir-Hinshelwood 'type of rate expressions 
which takes into account the effect of adsorption on rate of 
reactions. All three types of complex reactions, viz, indepen- 
dent, parallel and series reactions have been studied and the 
results are presented in the form of selectivity ratio as a 
function of modified Thiele parameter. Bor each type of reactions 
only one of the two reactions is considered reversible. The 
dimensionless parameters studied ares modified Thiele parameter, 
equilibrium constant for reversible reaction, K^, ratio of 
reaction rate constants, K^, K^p o , diffusivity ratios and 
adsorption equilibrium constant ratios with respect to reactant 

Diffusion with chemical reaction results in a set of 
non-linear differential equations of the boundary v al ue type. 

A new technique requiring much less computation time is dev el opr 
to obtain the missing initial conditions of the boundary value 
problem using the principle of invariant imbedding. 

The product distribution as indicated by selectivity 
ratio is markedly effected by diffusion limitations. The modif' 
Thiele parameter is one of the main variable influencing the 
value of selectivity' ratio. Bor values of ^ less than 0.2, the 
effect is found to be insignificant and at values of 0^ larger 
than five, the selectivity ratio attains a constant value. 

Except in two cases for independent reactions, selectivity ratio 
decreases with increase in 0^, The other variables which have 
a significant effect on selectivity ratio are K^, K^, 

Bg/D^ and K^p Q . Adsorption equilibrium constants and diffu- 
sivities of other reaction products are observed to have very 
small effect. 
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SIITOPSIS 


’'Effect of Diffusion and Adsorption on Selectivity of 
some Complex Reactions'” .A thesis submitted in partial fulfilment 
of the requirements for the Degree of Doctor of Philosophy by 
Avadh Behari Lai Agarwal to the Department of Chemical Engineering 
Indian Institute of ‘technology, Kanpur on October 1969. 

The effect of diffusion on the selectivity of some 
complex reactions in a cylindrical catalyst pore has been studied 
using Langmuir-Iiinshelwood type of rate expressions which take 
into account the effect of adsorption on rate of reactions. 

The following three types of complex reactions 
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are studied and the results are presented in the form of 
selectivity ratio, S^, as a function of modified Thiele 
parameter, 0^. The term selectivity ratio is defined as the 
ratio of selectivity with and without diffusion effects. 


f Ic^K^RT 


Modified Thiele parameter is defined as L)/ 


A v A^o' 

For each type of complex reactions only one of the two 
reactions is considered reversible. The dimensionless para- 
meters and the range of their values (given in bracket) studied 
are; modified Thiele parameter 0^(0, 1 - 10.0), reaction 
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equilibrium constant IC^(0.1 - 100.0), ratio of reaction rate 
constant for two reactions E^O.l - 10.0), K^p Q (0.01 - 100.0), 
ratio of dif fusivities for component B, 0, D and E with, respect 
to that of reactant A (0.25 - 4.0) and ratio of adsorption 
equilibrium constants for component B,C, I) and E with respect 
to reactant A (0.25 - 4.0). At a time only one variable is 
varied while keeping all other variables at a constant value 
of unity. 

Diffusion with chemical reaction give rise to a set of 
non-linear differential equations of the boundary value type. 

A new numerical technique is developed to obtain the missing 
initial conditions of the boundary value problem using the 
principle of invariant imbedding. 'The missing initial conditions 
are required in the present work as they are proportional to the 
overall reaction rates in the catalyst pore. This new technique 
is simple and takes only 3 seconds on IBM 7044 computer to 
obtain the two missing initial conditions for a given set of 
variables. Accurate estimation of the value of selectivity 
ratio by invariant imbedding technique is possible for values of 
modified Thiele parameter up to two, but for larger values, 
the results obtained by this technique are in error, because the 
increment used is not sufficiently small. Quasilinearization 
technique is used to calculate the values of selectivity ratio 
for larger values of up to ten. But this technique also 
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fails in many cases, especially for series, reactions, for ^ 
values larger than three. Computer time required by quasiline... 
zation technique to solve these equations for one set of 
variables varies from 15 to 60 seconds depending on 'the type of 
reactions and the values of the coefficients. 

The result shows that selectivity ratio of a complex 
reaction is significantly effected due to diffusion limitation 
in a catalyst pore. The magnitude depends upon the type of 
reaction, reaction mechanism and the physical and chemical 
constants influencing the system. The effect of different 
variables investigated are as follows; 

The selectivity ratio is found to decrease, except 
for two cases, with increase in modified Thiele parameter 
upto a value of five. Further . increase in does not 
effect selectivity ratio. Only for the values of Bg/D^ 
or Kg less than unity, the selectivity ratio is found to 
increase with increase in ^ in the case of independent 
reactions. For values of modified Thiele parameter less 
than 0.2, diffusion limitations are insignificant and 
selectivity ratio is close to unity. The selectivity ratio 
is found to increase with increase in the value of for 

independent and parallel reactions, while the effect is 
reversed for series reactions. YVhen the value of Kg is 
less than unity, the selectivity ratio is always greater 
than unity for independent reactions and the value increases 
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with, decrease in 11^ or increase in $2^, For such reactions ,wheir 

the value of 11 ^ is greater than unity, decreases rapidly with 

increase in the value of K^. For parallel and series reactions 

selectivity ratio decreases monotonicalljr as the value of K ^ is 

decreased. When Ej^P 0 is varied, the change in selectivity ratio 

is insignificant, if presented as a function of {2^. H 0 wever, 

when values are plotted as a function of ordinary Thiele 

paraaeter, it increases with increase in the K^p o value. The 

effect of increase in the value of diffusivity ratio, Dg/h^, 

is significant for all the reactions and the value of S R 

increases with increase in F-g/D^. Increase in diffusivity of 

component C increases selectivity ratio for independent and 

parallel reactions hut for series reactions, S R is found to 

decrease. The effect of changing the diffusivity ratio 

is insignificant for all the reactions. For independent reactions, 

the effect of change in the diffusivity ratio of two reactants, 

D„/D a , is similar to that of K , 

Jli A • 2 


The effect of change in the value of adsorption equilibrium 
constant ratios Kg/ K^, K^/l/ , K^/K^ and K-g/K^ on selectivity 
ratio is relatively insignificant. For independent and series 
reactions, selectivity -ratio is found to increase, while for 
parallel reactions it decreases with increase in the value of 
Kg/K A , K 0 A a and ICpA A * Comparatively, changes in K^/K^ &as 
maximum effect on selectivity ratio. For independent reactions, 
the selectivity ratio is found to increase with increase in the 
value of KgA A for % values up to two, hut for values of larger 
than, 2 . 3 , S R decreases with increase ip- the value of 
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C HAPTER - I 

I FZRODUCTIO I 

A. Background 

When tiie reaction takes place in a porous catalyst 
mass, the fluid reactant lias to diffuse inside the catalyst 
particle to make Hie effective use, of the surface area avai- 
lable in its pores since most of the surface area is inside 
the pores- as pore-walls. The reactants are chemisorbed and 
then react to form the chemisorbed product. The product mole- 
cules are then desorbed and diffuse out of the catalyst 
particle to the surface from which it is transferred to the 
bulk of the fluid. Each, of these steps offer some resistance 
and the rate of reaction taking place in a catalyst pore 
depends upon the magnitude of these resistances. In the past, 
expressions for the rate of reaction in a catalyst pore are 
obtained taking into account these resistances. If all these 
resistances are considered togather, the rate expression 
becomes too involved and in many oases it is not possible to 
obtain the rate expression in terms of measurable quantities, 
Bor this reason generally it is assumed that out of all these 
steps, one of then offer maximum resistance and accordingly 
all other resistances are neglected. Thus the rate of 
reaction equals : the rate of the slowest step. 

As a result of pore diffusion there exists a concen- 
tration gradient inside the pore and hence the reaction takes 
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place with differ.. nt rates at different positions in the 
catalyst pore, This effects the activity of the catalyst 
which is measured in terms of effectiveness factor which is 
defined as the actual rate of reaction in the pore divided 
by the rate of reaction if there is no concentration gradient 
due to pore diffusion. In complex reactions where more than 
one reaction is taking place, quite often, different reactions 
are effected in a different way and there is a change in 
selectivity, the term selectivity is generally defined as 
the rate of the desired reaction divided by the rate of 
disappearance of the key reactant. 

B. Literature Review 

The effect of diffusion in a catalyst pore was first 

treated by Damkohlar [l] , Thiele [2] and Zeldowitch [3]» All 

of them considered first order irreversible reaction and 

/ 

assumed that the lick s law of diffusion holds good. A new 
term, Thiele parameter 0 , was introduced which is qualitative!; 
the measure of the effect of pore diffusion on chemical 
reaction. The effectiveness factor was calculated as a 
function of Thiele parameter and it is observed that at 
higher values of Thiele parameter the effect is more while 
at values less than 0.5 the effect is negligible. Later the 
analysis was extended to include a wide range of reaction 
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order, pellet geometries and changes in number of molecules 
by Wheeler [4] and Weis and Prater [5]* 

In catalytic reactions, the rate expressions are not 
simple first or second order. To account for the effect of 
adsorption, Langrauir-IIinshelwood type of rate expressions 
are to be used. IP or the first time Chu and Hougen [6] 
considered the effect of adsorption on effectiveness factor. 
They analysed first order irreversible reaction A — j*B where 
the product adsorption is neglected. Roberts and 
Satterfield [7. 8] analysed the reaction, A + bB — » Products, 
assuming first and second order irreversible reaction kinetics. 
Schneider and Ilitschka [9, 10, 11 and 12] and ICao and 
Satterfield [13] considered first order reversible reaction, 
A^ B, and -showed that the retardation effect (d.ue to port,- 
diffusion) is more for reversible reaction than for irrever- 
sible reaction. 

Por complex reactions tile effect of pore diffusion 
on selectivity was studied by Wheeler [14] who designated 
three fundamentally different types of selectivity, since 
each is effected quite differently by the pore diffusion. 
Selectivity for independent reactions has been studied by 
Wheeler [14] for parallel reaction by Ostergaard [15] and 
Powlowski [16] and for series reactions by Carberry [17] 
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Weiz [5] and Wheeler [3.4]. All of them -have considered 
simple first or higher order reactions with simple rate 
expression. Recently Dwyer [13] lias studied selectivity 
effect by the competition among the three rate processes* 
that is, desorption, diffusion and surface reaction, in 
deuterium exchange of neopentane. 

G . The s is Obyect ves 

In the present work, effect of pore diffusion on 
selectivity of complex reactions is studied using Langmuir - 
Hinshelwood type of rate expressions which takes into accoun 
the effect of adsorption on the ra,te of reaction. Since 
there can be a wide variety of reactions, the following 
are considered for detailed study; 

1. Independent Reactions; 

A .0 + ■■ C 

E — D + C 

2. Parallel Reactions; 

2A -u==r 33 + G 

A D + C 

3. Series Reactions; 

A B + G 

B * D + C 

To consider the effect of reversibility, first of the two 
reactions in each case is taken as reversible. Por indepen' 3 
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and series reactions, single site mechanism is considered 
fox 'both the reactions where 3 and D react in the adsorbed 
state while C reacts in the fluid phase. For parallel 
reactions, first reaction is considered as reversible second 
order with dual site mechanism while for second reaction, 
single site mechanism is considered where D is reacting in the 
adsorbed state and G in the fluid phase. For all these 
reactions, rate expressions are considered where the surface 
reaction is controlling and reaction is considered in an 
infinite slab with cylindrical pores. Selectivity thus 
obtained is compared with the selectivity calculated in the 
absence of diffusion effects and the ratio of two selecti- 
vities, with and without diffusion, is presented as a 
function of fliiele parameter for various reaction parameters. 



CHAPTER - II 

MAh Ih i AT I C AL B0RMULA1I 0 I T CP T HE PROBLEM 
A* Independen t Rea ction s 

When two reactions with different react, ants &re 
talcing place on a catalyst, the reactions are called in- 
dependent reactions. Examples of this hind of reactions 
are hydrogenation of a mixture of aromatic and olefinic 
compounds, dehydrogenation of a mixture of napththenic 
and paraffinic compounds and cracking of a mixture of hydro 
carbons. The following reaction scheme is considered for 
detailed study. 

it 

B + 0 

( 1 ) 

I) + C 

If single site mechanism is assumed for both the reactions, 
where component 0 is reacting in the gas phase and surface 
reaction controls the rate of reaction, the following rate 
expressions can be obtained fox the components A and E. 


A 


E 


"1 


it 
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k 2 X Sl P E 


" r E ~ 


(i +<s iPi ) 

i 


( 2 *) 


Where and r 1? _ are the rate of formation of component A 


ana 


are the reaction 


and E per unit pore volume k-^, k^ 
rate constants for the adsorbed species in the direction 
indicated in equation 1, is the adsorption equilibrium 
constant for component i, is the partial pressure of 
component i, il is the equilibrium constant for the first 
reaction and subscript i stands for all ’idle components 
present in the reaction mixture vis., A, B, C, D and E. 
The rate of reaction of other components can be obtained 
from simple stoichiomet-.ry . 


r,, =~r. 
-u . A 


r D ~" r E 
r C = r B + r 3) 


\ J&J 


(?*>) 

(3c> 


Diffusion inside the pore t 

At steady state, a mass balance of the component 
diffusing in and out from a small thin shell of thickness 
dz at a. distance z from the pore^mouth becomes f ^ 


F'gtr 

(Rate of diffusion in the direction of z at z = z) - (Rate 
of diffusion in the direction of s at z = z + dz) - Rate of 
reaction in the shell. 
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ac, 

■ s -V«iT> 


ac, 


9 

a^c 


or D 


A 


2 

a^c, 


az £ 


[ ' SD A ( tf'- + d y? 


A 


a a^ 


)] 


Sdz r A (4) 

(5a) 


Where S is the cross-sectional area of the pore, is the 
concentration of component A and is the diffusivity of 


component A. 

The above equation is valid if transport of material 
is only due to diffusion and bull: transport due to pressure 
gradient is neglected. However, whenever there is change in 
number of moles, pressure gradient will exist within the pore. 
Otani et. al. [19] have studied the effect of pressure gra- 
dients on the effectiveness of porous catalysts and observed 
that for small capillaries where ICaucLsen. diffusion predominates 


the influence is insignificant, for a typical bidisperse 
catalyst pellet, they concluded that the effect can be importan 
only if the ct.ange in moles exceeds two. 

In general,: 


d 2 C, 

\ —? 

dz^ 


= " r i 


(5b) 


Substitution of the values of r A and r-g from equation (5) 


in equation (3) a gives 


2 

d^C. 


^ ah 


'B 


d 2 0 


- D, 


A 


da * 1 
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Integrating twice and substituting the boundary conditions 

c i( at * “ 0) - °io 

and 

dC. 

si d* z = 1) = o 

with L as the length of the pore, the following equation is 
obtained 

°3 = ho + <°A0 - V < 6 > 

Equation 6 is valid when diffusion constants are independent 
of the position within the pore, that is z and concentration. 
This condition is satisfied only for small pores where 
Knudsen diffusion is predominant. If we assume ideal gas 
behaviour, then 0 = P/B.I and 


P B ~ P Bo + Bg ^ P Ao " P A^ 


Similarly, p D = p ]Do + -g^(p. Eo 


(7a) 

(7b) 
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1C 1 K A [P A“ I ^Do+B^ ^ p Ao“ P A^ P Co + 3u (p Ao” p A )+ 

; jd o 


D A d P A 


i^ p e<t p e^ 


(8a) 


D . D 

1+K A P A +1 % Cp Bo + ^ (p Ao" P A^ +]I C [p 0o + D^ (p Ao“ P A )+ 


D^ P Eo" P E^ +K D^ P I)o + D d ^ P Eo" P E^ +K E P E { 


Similarly, 


d e d^p E 


k 2 i C E P E _ . . 

S 1+K A P A +K B ^- p B o + TC^ P Ao“ P A^ ^ +K c t p C o + Ih ^ p Ao“ p A^ + 


(81) 


E g ^ p Bo” P E.^ + K D^- P Do + I> e ^ P Eo“ P E^ +K E P E ^ 


Defining the following dimensionless parameters 

cp = p/p , rj = z/L, and K 1 = E/p q (where p Q is the pressure 

at the pore mouth) the equation 8 "becomes 

0?i>A- (a - /• V ( P- % V n- 


(uS+y'q). + 6 7 co ) 


(9a) 




Where ~ 


K 2 D E ( «/+ Y /( ? a + 6 9 b ) 


■ J ^K.Rf 

i of _ t# - A A — 


(9b) 


■ 


11 


1 


a 


9bo + Dg Y Ao 

' ^A „ , 2b 

^Co + D c ’ ^Ao + B r 9 Eo 


A 


ca , 


G 


« /= 1 + W a r; + p r; T k - -no ^ 


E 


S. + fi ^ (9 Do + £ 9 Eo ^ 


B 


( 101 ) 

(10c) 

(10 d) 


Y. 


Y 


K A p o (1 " 


K A p o ( X 

k l k a 


E 

A 


3 1 

D-m 


'A % 


K n = 


1 - kf Vo 


K 0 = 


k l K A 
k 2 IC E 


% % 

K 1 V 


% 

.*1 V 


(lOe) 

(lOf) 
(lOg). 
(10 1 '■ 


lie 'numerator and denominator 
divided ty 1 + K^p o > we ge u 


in equations 9a and 91 wlien 


B, 


d 2 cp. 


dr)‘ 


$n K- (a ~ ^ ? a 


D a ^ A 

4 CD ) (p - fY 9 


B c Y A 




(W + Y9 A + 6c P e ) 


r E y 


(11a) 


d2cp B _ & 1 

_ pr M 


B 


A 


9i? 


E 


dry 


k 2 d e 


(U + y9a + ^9^) 


(111) 


Here 


0m = ST ^Aa- 


El 


b a C1+V' 


(12a) 
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T 


__ 1 ^ 
d+J l/L p o 


-A P o 

1+r A p o 


ir 


K.p E, ' E, 

+ t|rV [ct f^r 

1 +K A P c li k . 1 A 


K-n 

+ |p(9--+ 


d e 

Lo t 9 Eo^ 


(1 - 


% ^A 

\ % 


5 d a , 

*P% 


, _ "A P o / hz 

5 = mt?,' 1 h 


K 0 D 
r A \ 


E 


% ?E 

-A % 


) 


(12b) 

(12c) 

(12d) 


with other constants remaining unchanged. 

Equations 11a and lib are to bo solved with the following 
boundary conditions 

9, (at '0 = 0) = 9, 

JO , ' 


r Ao 


cp E (at '0 = 0) = 9i 
(at T] = 1) 


(13) 


d°P A 




(at 'o = 1 ) = 0 


r Eo 

■ 19] 

d.i •' ■" = W 

ri:,e last boundary condition is obtained because no 

reaction and hence no diffusion can take place at the end 

of the pore. Recognising that the overall reaction xate 

within the pore equals the rate of mass transfer across the 

pore mouth, the rates of reactions of A and E will respectively 

l<h • _ dc? E 


be proportional to D ^ prp ^_q cml a.p 9=0 
selectivity with diffusion ex foe os 


and, therefore 


‘ rate of reaction of A 
S D = ra^*orYe¥c troll r>f B within the pore 
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The rate of reaction without diffusion effect can be ex- 
pressed as 

®Bo 9 Co 


"1 ^Apo ^°Ao 


) 


Ta = '^ Po ^' Ic i <£T 


(15a) 


1- 


ri2.r® f ojjo 

(i + p 


r B ~ 


^ E.cp. ) 
o ~*-£ iao ; 


and selectivity without diffusion 


< _ 

J 0 l'4n 


■\r t ( CD - ) 

.l A lAo. ' s 1 : : } _ 

lz 2 %E?8o 


(15b) 


(16) 


and selectivity ratio S R is defined as 


S T 


H ~ ^l/ S o 


Speci a l C a se ; 

Initially when no products are present, and also the 
two reactants arc present in equal proportion, then 

(17a) 


^Bo ~ 9 Co 


9-n = 0 

‘Bo 


p = co = 0,5 
r Ao 'Bo 


a =0.5 


J k 


Br 


(17b) 

(17c) 


p = 0.5 ( 


JO 

D 

1 


A + % ) 
B, J 


0 


0 


(17 d) 
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w = i_ + 5 &®g_ 


+ .5(24 

K A . D C 


E C , R K Dl \ 
r + ‘ 5 F T~i C 1 ?©) 
A V 


and Sq 





Therefore 



T) ,,.,'A ' 

A Jdr) rj=0 

V £% . " 

drf ' i—O 


(I7f) 

(I7g) 


with other terms unchanged. 

dcp, • 

To get selectivity ratio, the values of „ and 

d( p E • dr] r]=0 

-gpfp T) _0 are Jco "b e obtained. These are actually the missing 
initial conditions of the nonlinear simultaneous differenti al 
equations (equations 11a and lib) with boundary conditions 
^iven by equation 13. 


B . Pi a ral 1 e 1 R eaction s 

"Then two or more reactions with same reactants are 
taking place on a catalyst the reactions are called parallel 
reactions. Examples of this hind of reactions are dehydration 
of alcohols to olefins and ethers, decomposition of formic 
acid and catalytic cracking. The following reaction scheme 
is considered for detailed study. 
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If dual site mechanism is assumed for the first 
reaction and single site fcx the second reaction and also 
if it is assuried that the surface reaction controls the rate 


of reaction, the following rate expressions can be obtained 
for the component 3 and D 


1- r-* r - 0 2 „ jBjCs 
(i + i.K i p i ) 2 


(19 a) 


: 2 \ p A 


(1 + SXP,) 


(19b) 


Where r-g and r^ are the rate of formation of component B 
and D per unit volume and subscript i stands for all the 
components present in the reaction mixture. The rate of 
reaction of other components can be obtained from simple 


stoichiometry 


2 r B + r D 
r B + r 3 


(20a) 

(20b) 


Substituting equation 5 in equation 20 and integrating the 
equations thus obtained and substituting the proper boundary- 
conditions, the following equations are obtained. 


P A ~ p Ao ", ‘T)^ P 3 " P Bo ; “ t) A kp B" P Do 
' 3-q Ifv 

n = n -I — fn_— ' n_ ) 4- •a.*-- (ru— D— 


(p-n^l-no^ 


in ^ p D~ p Do^ 


PC = P Go + IT ^ 


(21a) 


(21b) 
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Tliese values of partial pressure when substituted in equation 
19a and eliminating r A from equation 51 the following 
equation is obtained 

2IL, Dj. 2Do ■ Dtj n2 
k l K A^ p io + D~ p Bo + D7 P llo - S: P B~ \ S I> 


Dg' a 2 P B 

.1 ~ ““**■ 
El 


O-Z 


X 

O p Co D n P 

J- ■ y 

* 23 ) 


1. 


'D, 


A 

X. 


1. 


D 


kO P 0o" B> - + Pti) ] 


S“ p Do + p; P B + D 0 y T>‘ 


B 


(22a) 


J-'-n 

{l +E A[PA 0 -r“( p B- p B o ^ P B~ p Bo ^ +P B P B + 


B 


V p Co" B^Bo" ^ P Bo + *4 Pb+ P C Pl)] + 


B 




D-i 


' K D P D j 


Similarly 


B d d 2 p 

Si 


B _ 


ZOb 

V 


D. 


h 2 I[ A^ p Ao + B~^* Pti - + tT p B" B. p B ) 


■■•2Db 


B 


35. 


P Bo + B^ P Bo~ D/B“ B/B 


dZ2 ^A [P A0 - 1^ (P B" S B0 ) - D^ p I)- p D0^ +K B P B + 

■ . D B- B ’ 


3). 


(22b) 


V P Co“ Bq P Bo” B“ P Bo + % P B + \ P D ]+K B P d] 

Introducing dimensionless pressure and length cp and P respect- 
ively in equations 22 the following equations are obtained. 

o 2 


d 2 9 


B 


dp' 


j> 

d X __ _ ~2_ Sjj__ ~ A - - 

dp 2 ( 9 b + 6 i 


||*| + l| 'Pb + 6 1 til' 


D 2Bq D-n 

^ \ '4 ( 0=1 ~ ^ ^ 


(23a) 

(23b) 
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where" 


Dt) 

^Ao + “hj % 0 + IT 9 Do 
^Co " D~ 9 Bo ~ $ 1)0 

1 + W“1 + ir Pi) 


( 24 a) 

(24b) 

(24c) 




iCW ^ iS uuu ^ A p o , ana U+£ a p 0 ) from the denominator .of 
equations 23a and 23h respectively, the>f flowing equations 
are obtained, \ 


: 



with, other coastants remaining ■unchanged. 

Equations 23a and 23h are to he solved with the following 
boundary conditions 

9 B (E = 0) = cp Bo (27a) 

9 d (t1 = 0) = 9d 0 (27b) 


^ dc ?B , ‘ ‘ 

dr) ^ ~ drl - (‘1 - 1) - 0 


(27o) 


fhe rate of formation of B and D will he proportional to 

d % J. dc P D , 

Bb~-^ jq_Q. and B^ V] _q respectively, and therefore, 

selectivity with diffusion effects. 

Rate of form ati o n of B within the pore 

Cppearance of A within the pore 


D A dt) 


T )=0 


o IS + S d? B 

B^ dTff)=0 1)^ *" dr) 


(28) 


r )=0 


Ihe rate of reaction without diffusion effects can be expressed 

aS ( ,„2 


l B 


Vlpo iS o ~ 

a 


„ ■ ¥A*i. 
r D - WTpyfr-^T 


(29 a) 


(29b) 


and selectivity without diffusion 
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K 


i-Vo^Ao - + p o ^ Vlo 5 


: t 1 + w )/ci + p 6 2 i K iW5 

and selectivity ratio will "be 

■ s D 

S H “ ■§“ 

0 

Special Case : 

Initially when no products are present 

'Pbo = fflo '= Tdo = 0 
cp Ao ss l (pure reactant) 

«i " 1 


h 


= o 


to. 


and S, 


Therefore S 


E 


K 2 K A P o 

i+z.f 0 (±~z;: 7 ) 

■ 1 +K a P 0 ( 1 + 2 K 2 ) 


d d 

®A 


■asr 




T)aaO 

^ * 4 * 

rt=0 D a Ar\ 


( 30 ) 


(31a) 

(31b) 

(31o) 

(31d) 

(31.) 

(31*) 


(31g) 


i)=0 


with other terns unchanged. 

To get selectivity ratio, the values of 
are to he obtained. 


n«o 


and 




T]s =0 
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C . Seri es Rea cti ons 

When the product of a reaction further reacts to 
form a new product ( s) * the reactions are called series 
reactions. Examples of this kind of reactions are dehy- 
dration of butanedial to form "butadiene, and dehydroge- 
nation of cyclohexane to cyclohexene and then to "benzene. 
The following reaction scheme is considered for detailed 
study. 

h 

A ' B + 0 

1 ; 1 («) 

kf 

B ► D + C 


If single site laeehrniism is assumed for 'both the reactions, 
where C is reacting in the gas phase and surface reaction 
controls the rate of reaction, the following rate expressions 
can be obtained for the component B and D 
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o tier components c an lie obtained from simple stoichiometry, 

" r A = r 3 + ( 34a) 

r C = * r A + r D (34b) 

Substituting equation 5 in equation 34 and integrating 
the equations thus obtained and substituting the proper 


boundary conditions, the following equations are obtained 

~ P^ ^ p B“ p Bo' ) ~ P^ p p“ p Po^ 

D-d 2P n 

+ l£ (p B*%o> + ^^ p B- p Po ) (35b) 


P A = P Ao 


P C “ p Co 


These values of partial pressure when substituted in 
equation 33a and eliminating r-g. from equation 5b the 

' ;V ' V 

following equation is obtained. 


P 


P T 


P T 


W P A0+ ]|W ft*T ^ P B" | P P 


% ^ 2 \ + jB + 

K ' lp Co~ P 0 ? 3o“ P G P Po P C P B + P Q P D' K 2 j 


Bb d-pg _ 

si ,2 p; “p^ Er 

[ 1 +K A ( p Ao + p^ p Bo + P^ p Po" P^ P B~ P A P P^ +il B P B 


V p c 


Pb 2P d P-g _ ^I)_ 


2D, 


o" P n p Bo 


p“ p po + p„ p b + p; p d' +ic p p p^ 




0 
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Similarly 


Dd d 


k 2 IC B P 3 

[ + ^(Py l0 + P]3 0 + p^Pdo" f~ P Pi “ Pj))+K-gPs+ 


t_ , % 23 D 2D T) 

V 0 (P Co- Dq P Bo~ DJT P Do + SJPb 4 ' p^ p p) +K D p j)] 


(36b) 


Introducing dimensionless pressure and length 9 and t] 
respectively in ecu '.rion 36, the following equations are 


obtained* 


r / jj -^tw P*d D*o 2Dj, 

[(a 2~ b a ^~ - r^z+l'V- irV-S 
v ;^r6r —rh-r k — 0 *- < **) 


C w 2 + ■ r 2?B + 6 2 >p B ) 




Talcing out l + K A p o from the denominator of equation 37a and 
37b and rearranging, the following equations are obtained. 


' 


la A . ^ 1 % D D “ 9 b A, 2D^ 

: f \ V ijV syV] 

^ ^ W 2 + ^2 % + * 29 !)) 


dr] : 


d*9 


D 


dp 2 

where 


5=^^ . ? B 

M ^2 ^ ru^Pf* ‘^r^gr 


= l^pT + 


K 


v k ■•= * 4 **> 


E 


1 K A p o 

Y 2 1+K a P 0 - a 


'D 


D 


2IL 


6. = 


K a P 


JA* 0 

1+ ^I P o 


( 


K 


E 


JD . ^0 s 

h + *T TT 5 


d a ' ~ T i D fl 5 


(39 a) 

(39 b) 

(40a) 

(40b) 

(40c) 


wi bh other constants remaining unchanged. 

Equations 39a and 39h are to he solved with the boundary 
conditions given by equation 27. The selectivity with diffusion 
effects is given b j equation 23, The rate of reaction without 
diffusion effects can be expressed as 

^Bo %o %, 


l B 


hK.p (cp. - 
1 A^o vy Ao 


“ jr 

% 


K, 


■) 


• (1 '+ p rt SL?,J 


(41a) 
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and selectivity without diffusion 


S 


0 



and selectivity ratio will be 

S; 




_D 

3 0 


Case ; 

Initially when no products are present 

V = ?Co = ?Do = 0 

1 (pure reactant) 


r Ao 


and 


a 2 

P 2 

bo g 

s o 


= i 
= o 

= l 
= i 


B, 


dqa 


B 


dr) 


T )=0 


Therefore S-o = -yr — 777- 

\ ~dp 


1. 


D 


dcp T 


b=0 + 'D a dr. H il=0 


with other terras unchanged 




To get selectivity ratio, the value of 
are to be obtai. 


T )=0 


. d . 9 . b ! 

dr) 


T]=Q 





(43a) 

(43b) 

(43c) 

(43d) 

(43e) 

(43f) 


(43g) 


and 
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C HAPTER - I II 

SOhUf IO II OF A SY SI 1 ELI 01? ITOI I -LIKEAR 
Ij OUITDARY VALUE BROBhEHS 


A . Int r o due t io n 


In a nonlinear differential equation when some 
conditions are given at one point while the remaining conditions 


at some other point, the equations are called the nonlinear 
■boundary value problem. These problems are much more diffi- 
cult to handle loth theoretically . and computationally, and 
there is no general proof of the exi stance and the uniqueness 
of the solution to problems of this type-. If all the condi- 
tions are known at one point, Bunge— Kutta or predictor- 
corrector method [20] can be used to solve a system of non- 
linear equations. Accordingly, one nietuod of solving non- 
linear boundary value problem will be to assume the missing 
initial conditions and solve these equations by one of the 
above mentioned methods. The final values are compared with 
the given boundary conditions , and the difference between the 
computed value and the true value is used to predict the 
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same as the true value, Shis difficulty becomes more severe, 

* 

if the equations arc of higher order or a large system of 
first order equations are to be solved. 

Several techniques have since been developed to solve 
linear boundary value problems of lirsr or higher order. 

!Ehes e are discussed by Fox [21] and Keller [22]. One of the 
methods is to convert these equations into a set of difference 
equations and solve them algebraically, bellman applied 
Newt on-B.aph.son s method to convert . noil linear equations into 
linear equations which can then be solved by the method of 
superposition or by finite difference method. The technique 
of converting nonlinear equations into linear equations is 
named as quasilinoarization and is discussed in his booh [23]. 
Quasilinearization technique has also been applied to the 
solution of a number of chemical engineering problems by 

Lee [24, 25, 26]. 

Quasilinearization technique nas been applied uo solve 
the nonlinear differential equations obtained in the previous 
chapter for fhiele parameter greater than 2 and is discussed 
briefly in the next section. For values of fliiele parameter 
upto 2, a newer method based on the principle of invariant 
imbedding is developed and is discussed in the next chapter. 


fhis method if 


MW. » s 


on quasi- 
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B. Qua a i line a^igauxog; 

Consider a system- of two second order differential 


equations 


,i2 


^ = f(x,y) = ^ 


dt 

£~X = g(x f j) = G 

dt 2 

with boundary conditions 
x(Q)= x Q 
y(0)= y 0 


f (1) = &(1> - 0 


(la) 

(lb) 

(2a) 

(2b) 

(2c) 


; tion at nth iteration is known, the 


If the value ox the iu.ua - 
value at n+lth iteration can be found by the recurrence 

relation 


_ 2 "X -p , \ B 1? 

iihtl = f( Xn ,y n ) + (V+l'Vi 5 ^ + (y n+l‘V 1>J n 


dt* 


(3a) 


d 2 y 


"B G 


£±i = g(x lx >y n ) + ( x n+l"* X n^ Cx,' + ^n+l y n^ By 


dt* 




XI 

(3b) 


and can bo solved by the methoc 


These equations one linear a 
of finite difference as dieoueeed in the neut section. 
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C . fi nite Diff erenc e Method * 

The left hand, side of equation 3 can he written as 


a % + i ( V , = 

dt 2 


“(At) 


Pi^L^n+l^^c+l^ ~ 2x n+l^h^ + x n+l^" b h-l' ) ^ 

(4a) 


- ririlh = n— JWW- ^n+fV + WWi 

dt ( At) (4.) 

To main algebraic equation simpler y n+1 is equated 

to y n and x a+1 is equated to ^ in equations 3a and 3b 
respectively. Substituting equation 4 in equation 3 and 
writing A for At tlie following relations are obtained: 


^ WW+? + vxriftA "as ^ih + i>= f ri‘A 


1 F 

X 


(5a) 


1 


G* 

/ + > (JL+ ££) y -i ( t J + -Xrin+l ( \+l) =G n" y n ( V ~Sy 

r y n+l (t h-l } v ^x^n+l v h ^ ^ 


(5h) 

If these equations are expanded for values of fc=l to H-l, 
equation 5a can he. written as 

■ x n+1 c\) + vi'V = V A^+l (t o’ 


-££*r>+^ 1 ; 


iri +a 2 WV + A 2 WV = h 

■ ' V . 1 ri . . - .. 

■ 1/ fl 
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+ a 3 WV + ^ 2 “WV = V 3 


KA-^^-iVA-i) '+ ~x 


^ 2 x n+l^%^ “ V-l 



u ne b ouncia r j c o nca bions are 


suasr irate 


where matrix A 



with. 


a = - (-£- + * 11 

h 1 A 2 • 




dx 

l S' 


_ _ AaA- + ALA) + 

1 h2 a* A 2 


) , k = 1 to IT-2 

1 


and 


V E n 


: n^ ?> 


Si? "o 


X i 

A 2 


■ ■ , . N > F 

V k~ ^n ” ^n^k "d x 


k = 2 to SI-1 


Similar expressions for equation 51 are 

5 A + 1 = W 

where B = F h-^ 1/A 

1/A 2 ^2 l/&2 

1 /a 2 t >3 !/a 2 


l/A 2 #l 3 I-l- 


^n+l 


T'il + l (X 2 > 

• '•# 

.. W 

* . ' A 

* 11+1 ( WJ 


( 12 ) 


(13) 


(14) 


(15) 


( 16 ). 
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with b., 


^N-l 


(‘ 


A 2 


O 

A. 


A 2 


+ 

•+ 


. OG n 

^ y ; 
de.£) 
?s y' 


j lv — 1 go IT** 2 


+ 


1 

A2 


and. 




- y (t-,) 

n - 1 - 


1G 

sy 







"AG 

Ty 


k = 2 to 1-1 


(18) 


(19) 


The method of solution consists of choosing a reasonable 
initial approximation of the function x(t) and y(t)j call them 
x (t) and y Q (t), ihe better values x-^(t) and y-^(t) are 
obtained by solving equations 8 and 14 » The improved value 
x-^(t) and y^(t) now becomes the initial approximation and 
the better values are obtained by again solving equations 8 
and 14. This process is continued till the two iteration 
gives nearly sue values for x(t) and y(t) . Equations 8 and 
14 can be solved easily by the : Thomas method [27] for tri- 
diagonal matrices. 
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CHAPTER - IV 
III VARIANT IMBEDDI NG 

S OLU TION OP A SYS TEM OF NONLIN EAR BOUNDARY VALUE 

PRO BLEMS 

A. Intro duo tion 

The principle of invariant imbedding was first used 
by Ambarzumian [26] and was later applied by Chandrasekhar [29] 
in treating the problem in radiative transfer. Bellman, 

Kalaba and Wing used them in the study of transport theory 
[30, 31, 32], A more detailed bibliography is given in a 
recent book by Lee [33] on Quasilinearization and Invariant 
Imbedding. Recently Srinivasan and ICoteswara Rao applied 
this principle to age dependent birth and death processes [34]. 
Bellman et.al. [35, 36] used this principle to solve nonlinear 
boundary value problem. Koenig [37] and more recently Lee [ 38 ] 
has demonstrated the usefulness of this approach for obtaining 
the missing initial condition of a second order differential 
equation. The method discussed by Lee uses linear inter- 
polation of values and is applicable only ror a set of two 
first order equations* Eor a system, of larger number of 
equations, Lee lias suggested a method in which a set of partial 
differential equations are to be solved. 



B . 'l!he sis Ob j e c tr vjcjs 

In the present worlc a new method is developed to get 
-the Hissing initial conditions of a system of nonlinear 
differential equations of the boundary value type. The 
method Is direct and requires significantly lees computation 

time . 

C . Mathema t ical 

Consider a system of nonlinear differential equation 

4Xj = f . (on , y l* J 2’ yn ’ ^ 


dt 

dt 


= g 1 (x 1 ,x 2 .x n , yq. y 2 - y n> t} 


i = 1 1 2 > D - 


( 1 ) 


with "boundary conditions 


Xi (0) 


c io 


( 2 ) 


y ± (t f ) = o 


with 0«tit f , t f is the final value of t. 

Do illustrate the approach, consider a system of 
only four non-linear two point boundary value equations. 
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g 2 ( x i * 


Jv Jo’ ^ 


with boundary conditions 
x 1 (0) = c 0 
X 2 (0) = d Q 
yi(%)= 0 

y 2 (-t f )= 0 

with 


lor solving e 4 uaUon S , consider a more general 


boundary conditions 

x^(a) = 0 . 

x 2 (a) = <?- (5) 

y^V = °. : ; 

y 2 (t f ) = 0 

Here a is starting value of tie xndependent 

variable t whloli varies from 0 to t £ . If ^(tj) and x 2 t 

d.u ig n0 longer a boundary value problem an 
are laaown- then r u is no ° . 

. . , . , lue problem where final conditions 

it becomes an xinual value pro 

• • -initial conditions y-,( a) and y 2 ^ a J 

are given. fhe missing ininal con 1 

- +0 ^ five direct method using 

wlll toe ol tailed * -» ^ ^ the missing 

■ invariant imbedding concept. It xs ap,are. 

-„ 4 , (a) for the system represente J 
conditions y 1 (a) a.n j z \*j 


x 
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equations 3 and 5 are not only tlie function of a, the 
starting value of tlie independent variable t, but also the 
function of initial conditions c and d, thus 


y-q(a) = r(c,d,a) = rtx-^a), x 2 (a), a] 
y 2 (a) = s(c,d,a) = slx^ia.) , x 2 (a), a] 


( 6 ) 


For a process starting at a + A , the missing initial condition 
y-^(a + A) can be related to y-^(a) by the use of Saylors 
series. 


y, (a + A) = y n (a) + y-,(a) A + higher order terms 


Where yj_(a) - ( ) t= . 


The value of y-^(a) caL1 " be obtained from equation 3 

y 1 (a) ~ gq[x-^(a) ? ? yq(a) ? y 2 ^®^ * a -l 

= g 1 [c,d,r(c,d,a) , s(c ,d,a) , a] (3) 


Substituting equations 6 and 8 in equation 7 and neglecting 
higher order torus, the following relation is obtained 

y 1 (a +A)=- r(c,d,a) + ^Cc,d,r;(q>d',a) , s(c,d,a),a]A 

(9a) 



d,a), s(c,d,a),a]A 
(9b) 

(c,d,a) ,a] A (9o) 
(c,d,a) ,a] A(9d) 


+ f [o,d,r(o,a, a ), 
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-Cte values of y-^a + a) and y 2 (a + A) can also be 

obtained from equation 6, thus 

y (a + A) = r [x-^a + A) , * 2 (a +A),a + A ] (i 0 ) 

y 2 (a +A) = s [ Xl (a +A), x 2 (a + A) ,a + A ] 

Substituting equations 9c and 9d in equation 10. 

yi (a + A) = r[c + ^[c^.r (c,d,a), s(c,d,a),a]A , 

d + f 2 [c,d,r(c,d,a), s(c,d,a),a] A, a+ A] 

(11a) 

yo ( a+ a) - e[o + fqLc ? d,r(c f d,a) , s(c ,d,a) ,a] A, 

d + fn[ o,d,r(c,d,a),s(c,d,a),a]A, a+ A ] 

2 (11b) 

The desired relationship can be obtained by equating 
equations 9a and 11a. 

1 ( 0 , a, a) + g 1 [c,d,r(o,d,a) ,s(o,a,a) ,a] A 

r [c+f, [=,d,l(o,a,a) a). ,a] A ,4rf 2 [o , 4,r (= ,4.a) , 

1 _ « 1 (12a) 

s( c , d, a) f a] A , a + A J 

Similarly ly equating equations 91 ®d 111 

s (c,d,a) + g 2 [c,a,r(o,d,a),s(c,a,a),a] A - 

[o,£,r(c,4,a) ,s(o,4,a) ,a]4 , d+f 2 [c , a, r(o ,d,a) , 

x ^ (12b) 

s(c , d, a) ,a] A ,a + A j 

- *■ „ q 12a and 121 are the required difference equations 

Equations 12a aau 

rvivt-qin the Hissing initial 
and can be solved airec « y u ’ 

conditions r(c,d,a) and s(c,d,a). 
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D. I Hus trative Example 

Consider tue following non-linear two point boundary 
value problem. 

d 2 x 1 (l-x 1 -x 2 ) 

dt 2 l+Kx^Xg) 

( 13 ) 

d 2 x 2 ■ -- (l-x 1 -x 2 ) 2 -x 2 
dt 2 [l+i(x 1 +x 2 ) ] 2 



( 14 ) 



given above is obtained when we consb 


der diffusion and reaction in a cylindrical pore for a parallel 


langmuir-HinshelYifood kinetic model 


reaction followir 


to obtain 


Replace 
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dy, 

dT 


(l-x 1 -x 2 ) -x 


= d 2 


(15) 


With more general boundary conditions 
x-^(a) = o 

x 2 (a) = d 

y^i) = o 
y 2 (i) = o 

Comparing equation 15 with equation 3. 
f 1 ( Xl ,x 2 ,y 1 ,y 2 ,t) = yi 

f 2 (x 1 ,x 2 ,y 1 ,y 2 ,t) =y 2 

g 1 (x 1 ,x 2 ,y 1 ,y 2 ,t) = G 1 (x 1 ,x 2 ) 

§2 ^ X 1 ,X 2 , *^'l * ^2. f = ^2^ X 1 ,X 2^ 


(16) 


(17) 


(18) 


Substituting equation 17 in equations 12a and 12b, equation 
18 is obtained. 

r(c,d,a) + G 1 (c,d)A =r[c+r(c,d,a) A ,d+s(c,d,a)A , a+ A] 

s(c,d,a) + G 2 (c,d)A =s[c+r(c,d,a)A f d-irs(c,d,a)A , a+A] 

If A is small, the following approximations can be made 
r[c+r(c,d,a)A , d+s(c,d,a)A ,a+A] = 

r[c+r(c,d,a+ A)A , d+s(c, d,a+ A): A ? a+ A] 


and 


■ 

ri j . a ( r> . fU sA A . a+ A "l == 


( 19 ) 


s[c+r(c t d # a)A ,a+s(c,d,a)A ,a+A] = 

s[o+r(c,a,a+A)' . 


A , a+ A] 
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and equation 18 can be written as 

r(c,d,a) =-G- L (e ? d)A +r[c+r(c,d,a+A) A ? d+s(c,d,a+A ) A ,a+ A] 
s(c , d,a) =-G 2 (c,'d) A +s[c+r(c , d,a+ A) A , d+s(c,d,a+A) ,a+A] 

( 20 ) 


With boundary conditions at t = 1 

r(c,d,l) = yqd) = 0 (21) 

s(c ,d,l) = y 2 (l) = o 

Divide the duration of t, 0 ^.t 1 in II parts with increment 

sS\ s o chat U ^ = 1 # 

To get r(c,d,l -A) and s(c,d,l-^), substitute 1-A for a 
in equation 20 thus, 

r(c,d,l~A) =-0 1 ( c fd)A ^ 22 ^ 

s(c,d,l-A) =-G 2 (c,d)A 

Put a = 1 — 2 A in equation 20 and substitute equation 22 to 
g et, 

r(c,d,l-2A) =^G 1 (c ? d)A +r[c-G 1 (c,d)A 2 ,d-G- 2 (c,d) A 2 , 1-A] 

s(c,d,l-2A) =-<} (c,d)A+s[c-G 1 (c ? d) A 2 ,d-G 2 (c,d) A 2 , 1-A] 

(25) 


( 22 ) 


let c-^ = c-G 1 ( c,d) A 
d-^ = d— G 2 (c,d) A 


(24) 






and since (from equation 22 by replacing o and d with c- L and 
respectively) 

r(c 1 fd 1 ,l- A ) = ]* G i( c i>di)A (25) 

( JT% - A .1-A } 22 — G^ ( C-1 • d-1 ) A 
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equation 23 simplifies to 

r(c,d,l-2 A) =— [ G-^ ( c , d) + G-| (c^,d-^) ] a. 
s(cjdyl— 2Zi) =--[&2(o*(i) + GgC ] A 


( 26 ) 


Similarly when, a = 1-3A 

r(c,d,l~3A>) =-G 1 (c, d)A + r[c+r(c,d,l-2fi) A , 
d+s(e,d,l-2A)A , 1-2 4] 
sfc,d,l~3A) t d)A + s[o+r(c,d,l-2A)A , 

d+s(e,d,l~2A) A , 1-2A] (27) 

Substituting equation 26 for r(o ,d,l-2A ) and s(c,d,l-2A) 


and defining 

c 2 = o-[& 1 (c,d) + & 1 (c 1 ,d 1 )] A 2 
d 2 = d-[G 2 (c,d) + G 2 (V d i)] 


(28) 


we get 

x(c,d,l~3A) = - G n (c,d)A + r(c ?) h,l-2A) 

x -- ( 29 ) 

s(c,d,l-3 A.) = *’* Ggfcjd)^ + s( Cg? dg » 1-2 A) 

The last terms of the above equations can be calculated from 
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r(c,d,l-3 A) 
s( C y d y 1-3 4i.) 


[& 1 (c jl) + & !C c i» d i^ + G l^°2» 

[G (c ,d) + & 2 (o x , d 1 ) f & 2 ^ c 2 ? d 2^.^ 


(31) 


The procedure when extended upto a = 1 - H 1S > 

_ o) , the desired missing initial condi'cions are obtained. 


y x (0) = r(c 0 ,d 0 ,0) = - 

y 2 (o) = s ^ c o ,d o , °^ 



(32) 


‘ ' QT , fl a _ a. The values of c^e,. . . 

since at a = Oj c - c Q and a o 

, , i d 0 xe calculated from the following relatione 

and d 2 ,d^.. axu 

, i u r/a A form of equations 24 and 28* 
which are the generalized ioim ox e<i 


X-l 

C T =C -[6,(0, 4) + X) G l (o i> d i^ A 
1 1=1 

1-1 

d I =d -[& 2 (c,d> + & 2 (c i? d i )]^ 3 


(33) 



from 2 to H- Hie values ox c j _ 

L successively as explained in the 

rever it is to bo noted that a new 
xlatad from the last computed value 
n ooMirsfit of the solution 


value of c 
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decrease in the value of increment,, but. ire amount of compu- 
tation and therefore computation time required for solution 
is also increased substantially* It can easily be shown 
that 2 x( 2) ” is the total number of c and d values calculated 
to get the appropriate values of c J s and d ? s for substitution 
in equation (32), but the solution is explicit and inter- 
polation or trial and error are not involved. 

The computational procedure to get c^ is explained 
in fable I, and the procedure is continued upto c^. The 
number in parenthesis in the fourth column indicates the 
step in which a particular value is obtained for substitution. 
By comparing equations 32 and 33 it can be shown that the 
missing initial conditions are 

y 1 (°) = ( ) / & 

y g (0) = (djy d Q )/^ 
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l 1 ABLE I 

Computational .Procedure to Obtain the Value of 


Step 

No, 

Quantity 

Calculated 

Equation 

. n.qe.l _ 

Values ^BsTi’fcutedT dumber 
__ in equation divisions.? 

(i) 

C 1 

24 

Oq for c 1 

(ii) 

. O' 

2 

28 

Cq for c and c-^(i) 2 

(ill) 

C 1 

24 

Cg(ii) for c 

(iv) 

°3 

33 ?i~3 

Cq for c,c-^(iii) and c 2 (ii) 3 

(v) 

.O-j , 

24 

c^(iv) for c 

(vi) 

C 2 

28 

c^(iv) for c and c-^v) 

(vii) 

°1 

24 

c 0 (vi) for c 

( viii) 

w /i 

H- 

33? i=4 

c^, for o? c-^vii) c^Cvi) 
and c^(iv) 4 

(ix) 

°1 

24 

c^(viii) for c 

(x) 

C 2 

28 

c^(viii) for c and c-^Cix) 

(xi) 

0 “7 

i- 

24 

c^Cx) for c 

(xii) 

0~ 

3 

33 ?i=3 

c^(viii) for c, c-^Cxi) 
and Cq(x) 

(xiii) 

°i 

24 

c^(xii) for c 

( xiv) 

C o 

c. . 

28 

c^(xii) for c and c-^(xiii) 

(xv) 

c i 

24 

c 2 (xiv) for c 

(xvi) 

c 5 

'35,1=5 

c 0 for c, (xv) , c 2 (xiv) , 

c^(xii) and c^(viii) 5 
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E . Results and Disc ussion 

The method presented in this paper to solve a system 
of non-linear boundary value equations using invariant imbedding 
concept is direct involving less computational error and time, 
fable II gives the values of y^(0) and y.^(0) obtained by 
solving equation 13 and 14 with c q = d^ = 0. 

TABLE II 

Missing Initial Conditions y-^(0) and y 2 ( 0) for CQ=dQ=0 


11 

A 

yq(0) 

y 2 (o) 

8 

1/8 

0.642778 

0.432815 

9 

1/9 

0.640918 

0.429448 

10 

1/10 

0.639431 

0.426793 

True 

value 

0.626181 

0.404287 


The computation was carried out on IBM-7044 and required 
approximately l|>- seconds for calculation with A=0.1 (R=10) . 
With each increase in the value of N by 1, the amount of time 
for computation doubles, thus for E=20, the time will be about 
-§• hour. To save computation time and still get a good accuracy 
a polynomial equation of the form given below may be used and 
the final values obtained by substituting A = 0. 
y x ( 0) = A ] _ + b iA + c 1 A 2 + 
y 2 (o) = a 2 + b 2 a+ c 2 + 


( 34 ) 
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A three constant polynomial with y^(0) and y^ (0) 
values calculated for A = 1/8, 1/9 and 1/10 gives the 
following values for the missing initial conditions. 

y 1 (0) = 0.626124 
y 2 (0) = 0.404393 

These values are very close to the true values which 
are obtained by the method of trial and error. 
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CHAPTER - Y 
RE SHIPS AIU) DISCITSSIOIT 

Results are ’ presented in the tabular form in 
Appendix A-l throu gh A- 29 and the plotted values are shown 
in Pig. 1 to 29. The following dimensionless variables are 
studied: 

- Modified Thiele parameter, 

- Reaction equilibrium constant, K-^ 

- Ratio of reaction velocity constants, 

- K A P a 

- Ratio of diffusivities, Bg/D^, D^/D^, , B^/D^ , and 

- Ratio of adsorption equilibrium constants, Kg/K^, K^/EC^, 

KrAa a2ld KgAk 

The number of variables studied, excluding modified Thiele 
parameter, are eleven for independent reactions and nine for 
parallel as well as series reactions. By combining the 
variables it was possible to reduce the total number of lumped 
variables to a minimum of six, but this method was not followed 
because of its inability to present the results in a clearcut 
way. On the other hand, the variation of these lumped variables 
(which are function of the variables studied) have little 
meaning as the variation in one will effect the value of idle 
other. Por these reasons, the variables are kept separate in 
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this study and only one variable is varied at a time. This 

is advantageous as it shows clearly the effect of the variation 

in one particular variable while all others are hept unchanged- 

However, p^, contains which is also considered as an 

independent variable, ihe inclusion of K^p o in- P^ is justified 

on the basis that it provides a better picture of the complex 

kinetic expressions and a very important generalization is 

possible, viz. for values of greater than five, in general, 

the selectivity effect becomes constant and further increase 

in does not change the selectivity of the reactions. Also 

selectivity effect for 0^. values of less than 0.2 is, in general, 

negligible. She use of lumped variables by many authors [6, 8 

and 9] to present activity of porous catalyst for single reaction 

is useful since it was possible to reduce the total number of 

lumped variables in those simple cases to one or two, in 

addition to modified Thiele parameter, and the presentation of 

the results was not difficult, A problem can be solved by 

finding the values of each of these lumped variables and the 

activity value is obtained from the proper curves. With six 

lumped variables, in addition to modified Thiele parameter, 

6 

at least 729, that is, 3 curves are required to present the 
effect on selectivity ratio as a function of pl» for only three 
selected values of each of these lumped variables for any 
given reaction, obviously a formidable task of doubtful merit. 



The effect of diffusion is presented in the form of a 
dimensionless selectivity ratio, S^, as a function of modified 
(Thiele parameter, , At a time only one variable is studied 
and the values of all other variables are Icept constant at unity « 
She range of the values of each of the variables studied is 
given in Table III. 


Ta ble III 

Range _of_ jdie^ Values of the Va ria b les S tudied 


Variable 

Range | 


0.1 - 10. 


h 

0.1 - 100. 


K 2 

0.1 - 10. 


Vo 

0.01 - 100. 


■Ohv W h/h 

0.25 - 4. 


KbA a , vh 

0.25 - 4. 


KcA a , Vh 

0.0 - 4. 

i 

, J 


for clarity of presentation true reaction rate constant 
and adsorption equilibrium constant are shown separately while 
defining the reaction rate constant, but the observed reaction 
rate constant is simply a product of the true surface reaction 
rate constant and the appropriate adsorption equilibrium 
constant. 
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/kjK.BI 

Modified Thiele parameter, is defined aa XV ^ y 4 qf “ p ‘ r ' y * 

A A o 

which is also a function of K^p . While studying the effect of 
variation in the value s of E^p o , the value of modified Thiele 
parameter also changes* In reporting the results, the variation 
in modified Thiele parameter due to change in the value of K^p 
is not considered. However if the net effect on the value of 
selectivity ratio due to change in the value of E^p o is required, 
the consequent change in (due to change in E^p Q ) is to he 
calculated and the selectivity ratio value is to be obtained at 

this new value of gL, In all such oases, the selectivity ratio 

11 

values are also reported as a function of ordinary Thiele 

* t— ■nm 

and is independent 
of K^Pq, The values of ^ as a function of 0 are given in Table IT. 

Very small, values of K.p mean that the modified Thiele 
parameter becomes nearly equal to ordinary Thiele parameter and 
the rate equations are nearly simple first or second order 
provided adsorption equilibrium constants fox other components 
are not large. When one of the product is not adsorbed at all, 
K./K. for that component is zero. Of special interest is the 

1 i. • ■ , 

case for component 0. When ’K^/K^—0 , the component C is not 
adsorbed at all and reacts in the gas phase only* When the 
value of reaction equilibrium constant K-^ is greater than tea 

the reverse reaction is negligible and hence the reaction can 

be treated as irreversible. 


par ame ter, p , whi cii is d e f ine das 1”V -- * * — 

^A 
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jPai;\r. _:vr 


Modifie'd Thiele Parameter , as a Eimc ti on of Ordinary 

Thiele Parameter, $, for Selected Values of K '.^p Q 



!\ K A p„ 

Aro 

V X 

0.01 

0.1 

1 . 

10. 

100. | 

0.1 

0.099 

0.095 

0.071 

0.030 

! 

o.oi i 

i 

1 . 

0.995 

0.954 

0.707 

0.302 

0.099 

. 10. 

1 

i— » - — 

9.95 

9.54 

7.07 

3.02 

0.995 


In all cases except two, the selectivity?' ratio decreases 
with the increase in modified Thiele parameter. Only when the 
values of D-^/D^ and of are less than unity, the selectivity 
ratio is found to increase with the increase in modified Thiele 
parameter. The change in selectivity ratio continues up to 
modified Thiele parameter, value of nearly five. For values 
°f S^ter five, the selectivity?' ra,tio becomes nearly 

constant, since at these values most of the reaction is over in 
a catalyst pore and further increase in modified Thiele parameter 


do not change the product concentration appreciably and therefore 
the selectivity ratio approaches a constant value. However, 
further increase in the values of will continue to decrease 
the activity of the catalyst. The concentration at the end of 


the pore has been calculated for parallel reactions for different 
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values of modified Thiele parameter keeping all other variables 
at unity and tlie results are presented in fable V. It is clear 
from tlie values given in Table V that the percent composition 
of reactants A and B for a value of greater than five is 
close "to zero, while that of the products C and D is nearly 
constant. 

Table V 


Mole Fractions o f .Di ffe r ent Componen t s a t Ihe P ore End, as a 
Function of Modif ied Thi ele Par a m e ter >. ..for* P arallel Reac ti ons 



1. I 0.441 0.060 0.280 0.219 

- * . . 


2. j 0.171 0.046 0.415 0.368 

3. 0.075 0.028 0.461 0.435 

4. 0.037 0.016 0.432 0.465 

5. 0.019 0.010 0.490 0.481 

7. 0.005 0.004 0.498 0.493 

10. I 0.001 0.001 0.500 0.498 



Method of invariant imbedding is used to calculate 
selectivity ratio lor different values of modified Thiele 
parameter up to a value two (in some oases up to a value of 


3 also). For values of 0^ greater than three, the values of 
selectivity ratio obtained by above technique is in error, 
because the increment used in computation (A - 0.1) is too 



large for these conditions. Quasilinearization technique is 
used to calculate the selectivity ratio for the values of 
modified Thiele parameter greater than two. This technique 
also does not converge in some cases for large values of 
9 ^( 9 ^ 10) in parallel reactions and for values of greater 

than three (and in some cases for = 3 also) in the case of 
series reactions. The time taken for obtaining one set of 
value s by invariant imbedding technique on IBM 7044 is less 
than three seconds. The time taken by quasilinearization 
technique depends upon the ‘type of reaction, the values of the 
constants, and the ac curacy desired. This time varies in these 
computations from 12 to 60 seconds, which is much larger when 
compared to 3 seconds required by invariant imbedding technique. 
The curves are extrapolated in the cases where the solution is 
not obtained and the extrapolated portion of the curves is 
shown by dotted lines. Mo value is reported in the tables for 
such cases. 

The discussion is based on reactions 
A ?=* B + C 
B > D + C 

and the selectivity ratio is defined as the ratio of selectivity 

with and without diffusion effects. The selectivity is defined 
as the ratio of the rates of disappearance of reactants A and B, 
that is, Bor the reaction system considered, the 
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selectivity without diffusion limitation is equal to 

For irreversible reaction, that is, K-^ greater than ten, 
if the dif fusivities of both the reactants A and. E are same and 
also, if both the reaction rate constants are equal in magnitude, 
then there is no effect of diffusion, because diffusion limi- 
tation will not change the ratio of reactant concentrations, aid. 
hence their relative rates* This is shown in Eig.l for = 100, 
where the value of selectivity ratio remains nearly unchanged. 

As the value of is lowered, the first reaction becomes more 
and more reversible and the net rate of the first reaction is 
reduced due to increase in the product concentration as a result 
of diffusion limitation. However the rate of irreversible 
reaction is not changed and therefore selectivity ratio is 
reduced. For highly reversible reactions, the reduction in 
selectivity ratio is quite significant [Appendix A-l] . 

Ihen is varied, the reaction velocity constant for 
the second reaction is changed. If second reaction is faster 
than the first and if the component E cannot diffuse equally 
faster, then the rate of the second reaction does not increase 
proportionately, since the surface concentration of reactant E 
keeps on decreasing due to diffusion limitation. Accordingly 
if the second reaction is slower than the first (Kg greater 
than unity) , eventhough selectivity is increased but the selecti- 
vity ratio is decreased with the increase in modified Thiele 
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parameter, 0^. Conversely, if the second reaction is faster 
than first (Kg less than -unity), the selectivity decreases 
hut selectivity ratio is increased -with the increase in & 
jFig.2 and -Appendix A- 2] . When Kg = 0,1, the value of. selectivity 
ratio in strong diffusion region (0^ = 10) is 2,43* If the 
reaction is considered to he irreversible (say K^ = 100), with 
Kg = 0.1 the calculated value of selectivity ratio is 3.04 
(not reported in tables). Furthermore, if irreversible reaction 
with pimple rate expressions (K^ = 100 and K^p Q = 0.01) is 
considered, the calculated value of selectivity ratio is 3.18 
(not reported in tables) . This value compares favourably with 
the value of 3.16, reported by Vfheeler [4] for such reactions. 

At higher values of K.p Q , the value of modified Thiele 


parameter decreases, which in turn results in a larger value 
of selectivity ratio. Fig. 3 is a plot of selectivity ratio 

versus modified Thiele parameter for different values of K.p 

A o 

without accounting for the consequent change in the value of 0^. 
The curves are so close that only two curves for the K.p values 

ii U 

of 0.01 and 100. are shown. When presented in this way, ten 

thousand fold increase in the value of K.p„ did not alter the 

A o 

value of selectivity ratio by more than two percent [Appendix A~3] 
The selectivity ratio values are slightly more for higher values 
of K a p q . The effect of varying can be seen more clearly, 

if (1 + k a p o ) is tahen out of the modified Thiele parameter and 
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•the values of selectivity ratio presented as a function of 

' fc,K A MI ' 

ordinary Thiele parameter, 0 - L V , for selected 

A 

values of K^p , keeping the value of all other variables at unity. 
When presented this way [Table VI] the spread between the curves 
is increased to a maximum of twenty percent and the selectivity- 
ratio increases 110:10 tonic ally with increase in the value of K.p 

■1 x O 

at any given value of 0. The increase in the value of K.p 

A O 

corresponds to increase in the value of denominator of rate 


equations, which results in the decrease in reaction rates. This 
is then analogous to lower effective values of Thiele parameter 
corresponding to higher selectivity ratio. 


T abl e VI 


Selectivity Ratio, S^, as a Function of Ordinary Thiele Parameter 



0t : 

:or Selected 

Values 

of *#0 


! K.p 

A 0 






X. 

i 

, : 0,01 

0.1 

1. 

10. 

100. 

0.1 

1.000 

1.000 

1.000 

1.000 

1.000 

1 . 

0.946 

0.955 

0.982 

0.999 

1.000 I 

2 . 

: 0.847 

0.855 

0.909 

0.990 

1.000 ! 

5. 

0.021 

0.824 

0.858 

0.968 

0.999 | 

4. 

0.016 

0.318 

0.85S 

0.959 

0.997 ! 

- • t. 

H 

O 

• 

ITS 

H 

CO 

* 

O 

0.816 

0.824 

0.842 

0.960 ! 

i 

" i 



59 


When dif fusivity of component B ox 0 is increased, the 
reaction product diffuses faster and therefore the concentration 
of product molecules are lowered. Consequently the net rate of 
first reaction is increased which results in an increase in 
selectivity ratio [ fig. 4 and 5 and Appendix A-4 and A-5]* 

The effect of changing the diffusivity of component D 
is negligible as a sixteen fold change in h^/D^ changes selecti- 
vity ratio only slightly. This is because the second reaction 
is irreversible and the concentration of component D has little 
effect on the reaction rates [fig. 6], Decrease in the diffusi- 
vity of component D increases its concentration,- thus lowering 
the rates of reaction due to the increase in "die value of 
denominator in rate equations* This effect is then analogous 
to the effect of decrease in the value of , that is, lower 
values of D-^/D^ results in higher values of selectivity ratio 
[Appendix A-6], 

The effect cf changing dif fusivi^ of reactant E with 
respect to reactant A [fig. 7 and Appendix A-7] is similar to 
that of Kg as shown in fig. 2* When the diffusivity of component 
E is large, the decrease in concentration along the pore will 
be small, and hence the rate of second reaction will increase 
resulting in the decrease in selectivity ratio. 

When the value of adsorption equilibrium constant for 
reaction products B, C ox D are increased, the selectivity 
ratio increases. Since increase in the value of adsorption 

I I 1 n . ' iff . Slill 











r 


liiBliiil 


1-6 


PialliiiJii 




1-2 



0-0 

TiiWi 


. ! “ t *' 1 i - 1 1 

"re ?■::< ' 40 io-o 

l/V / . 1 '■ : : 


ig7- Selectivity ratio, S r, as a function of modified 
Thiele parameter,(J)^, for seleqted values of Dr^D/ 
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equilibrium constants increases tlie value of the denominator 
which in turn decreases the rate of reactions and consequently 
increases the selectivity ratio [Pig- 3, 9, and 10] as explained 
earlier- All the curves lie in a very narrow band a n d the 
sixteen fold change in the value of adsorption equilibrium 
constants result in less than two percent change in selectivity 
ratio for component !3 and D and a maximum change of three percent 
for component 0 [Appendix A-8, A-9 and A~10]. 

The effect of changing the value oi adsorption equilibrium, 
constant for the reactant component E is sane as that of product 
component B, C and D for lower values of modified Thiele parameter 
that is, for higher values of the selectivity ratio is 

larger. However, the effect is reversed when the value of 
modified Phi el e parameter is greater than two [Pig. 11 and 
Appendix A-ll] . 


I'he computed values of the rate of disappearance of 
reactant A and E shows that the reduction in the rate of 
disappearance of reactant E is always greater than that of 
reactant A as the value of Kg/iC^', or increases 

from 0.25 to 4.0 for all values of 0^. However, tins trend 
is found to be valid while varying from 0.25 to 4*0 onlj 

upto a modified ihi cl e parameter value of less than two. Por 
higher- values of rhiele parameter with increase in tlie value 
of ICg/K^, the reduction in the rate of disappearance of 



0-1 0-2 0-4 1-0 2-0 4-0 10-0 


Fig.8-Selectivity ratio, Sr, as a function of modified 
Thiele pa r ameter .(p^ , for selected values of Kq/K; 









Fiq.11 -Selectivity ratio, Sr, as a function of modified 
Thiele parameter ,0^, for selected values of K F /K A 
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reactant A is greater titan that of reactant E, and this reversal 
is then responsible for the decrea.se in the selectivity ratio 
with increase in tie value of Kg/EC^ for 0L values greater than 
two. For large values of $L, the concentrc.ti on of reactant E 
is close to zero towards the end of the pore because second 
reaction is irreversible aid this is probably responsible 'chat 
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B . Par allel J&eacjy.oiis 

The discus sion is based on the reactions 


2A ^=-: B + C 
A“ — * D + C 

and the selectivity ratio is defined as the ratio of selectivity 
with and without diffusion effects. The selectivity is defined 
as the ratio of the rate of formation of product B to that of 
disappearance of reactant A, that is, r-^/(-r^). Por the 
reaction system considered, the selectivity without diffusion 
limitation is ecjual to I"*"2 x K A P o /(l-hK A P o (l+2K 2 )). 

The variables studied for parallel reactions are same 
as listed in Table III except that the component E Is not 
present. 

When the first reaction is highly reversible corresponding 
to low values of equilibrium constant li selectivity ratio is 
decreased due to rapid increase in 'hie concentration of products 
C and P as a result of diffusion limitation. The concentration 
of component b is first increased and then starts decreasing 
as more and more reactant A is converted to products G and D 


to maintain the equilibrium of first reaction. Thus for a given 


value of Thiele parameter as the value of K 2 is lowered the 

rate of formation of B decreases rapidly while that of D 
increases slightly due to diffusion limitations. This results 

in the decrease of selectivity with diffusion limitations, 


which results in 



Hues of selectivity ratio [Pig. 12 
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3?or a given value of the change in Kg implies the 
change in the value of rate constant for the second reaction 
without effecting the value of tire rate constant for first 
reaction. When the value of Kg is decreased, the rate, of second 
reaction is increased which decreases the reactant concentration 
and increases the product concentration much more, the net 
rate of first reaction is therefore reduced much more (since 
it is second order reversible reaction) than that of second 
reaction resulting in lower selectivity ratio. Conversely, if 
the value of Kg is increased, the second reaction becomes 
slower and the decrease in reactant concentration as well as 
increase in product concentration along the pore is reduced. 

She net rate of first reaction is therefore increased much more 


than that of second reaction resulting in higher (but always 
less than unity) selectivity ratio [Kig.13 and Appendix A-13] • 
When the value of Kg is large, the rate of second reaction is 
low and therefore even for large values of modified Thiele 
parameter, conversion of reactant A to products C and D is 
not complete and concentration profile changes with change in 
the value of $, T . As a result of this, the value of selectivity 
ratio does not become constant for p^ as large as ten, when 
Kg is greater than two. However when second reaction is very 
fast the conversion of reactant A to products C and 3) is 
achieved for values of as low as four (Kg = 0.1) and 

further increase in does not alter the value of S-, 

, 

Wmlmmm 


. 


R 
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K A P 0 effects tne rate equations in three w ay s .. This 

appears as a coefficient in tlie rate equation for first 

reaction, in nodified Thiele parameter and in the denominator 

of the rate equations. Pig. 14 is a plot of selectivity ratio 

versus modified Thiele parameter for different values of Ii.p 

A^o . 

without considering the effect of the variation of IC p on jX T . 
When presented this way ten thousand fold increase in the value 
of K a p q alters the value of selectivity ratio to a maximum of 
twenty percent at a ^ value of ten. furthermore, the value of 
selectivity ratio are slightly higher for low values of ^ and 
K a P 0 in comparasion to those at low values of ^ but large 
values of K^p Q , Ibr $2^ values greater than one, the selectivity 
ratio increases with the increase in the value of K.p at any 
given value of fL [Appendix A-14]. 

To show the effect of varying ihp more clearly, the 
values of selectivity ratio are presented in Table VII as a 
function of ordinary Thiele parameter, for selected values 
of K a P 0 » When presented this way, the spread between the 
curves is inc : eased quite considerably and the selectivity 
ratio is found to increase monotonically with increase in the 


value of K a P q at any given value of 

The selectivity ratio values for ordinary as well as 
modified Thiele parameter greater than three and for E a P q value 

of 0*1 are not conforming the general trend as can be seen from 
Table VII and Appendix A-14. No particular reason for this 
behaviour is evident. 




as a function of modified 
selected values of K^p t -,. 
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Table VII 

Sel ectively Patio , S R , as a function of Ordinary 
Thiele Parameter, p', for Selected Values of K.n 

* »• . _ A-^n 


j x 

! Gf l\ 

! 4 X 

t*-— — - - .... - v ■ ... ' , „ ■ 

0.01 

! 

1 

0.1 

1. 

10. 

100. I 

i 

0.1 

0.997 

0.997 

0.998 

1.000 

1.000 

j 1. 

t 

0.777 

0.788 

0 *864 

0.965 

0.996 

! 2. 

0.567 

0.591 

0.709 

0.897 

0.984 

j 3. 

0.515 

0.518 

0.618 

0.850 

0.966 

! 4. 

i 

0.498 

0.494 

0.575 

0.775 

0.945 

I 10. 

i — 

0.485 

— 

0.551 

0.617 

0.815 


Wlien dif fusivitie s of component 33 or 0 are increased 
the reaction product diffuses faster and therefore the con- 
centration of the product molecules is lowered. This results 
in lowering tne reaction rate of the reversible step and there- 
fore the net rate of first reaction is increased, resulting in 
an increase in selectivity ratio [Pig.15, 16 and Appendix A-15 
and A-16] . 

of component D 
the small difference 
n ihe value of the 





D b /D a =io' 


■ , fVJ , , . 

Selectivity ratio, Sg,os a function of modified 
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<p M 

Activity ratio, S Rj as a function of modified 
parameter ,(p M , for selected values of Ch/D A 


0-9 


0-8 

c r 

if) 

0-7 
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denominator of tiie rate equations as a result of tire change in 
the value of Dp/D^, The effect is similar to that reported for 
independent reactions. 

The effect of increasing the value of adsorption equilibrium 
constant for reaction products B, C or D is to increase the 
value of denominator of the rate equations. The reaction rates 
are therefore reduced, hut the reduction for first reaction is 
much more in comparision to second reaction, since it is a second 
order reaction. Consequently the value of selectivity ratio is 
lowered. The spread in the curves is maximum when Kg is varied 
and is minimum for the variation in Kg, The component C parti- 
cipate in the reverse reaction and its concentration is equal 
to the sum of the concentration of component B and I) and there- 
fore the effect of varying Kg is more than for varying Kg or Kg. 
Even though component B also participates in reverse reaction, 
hut its concentration is always small [Table V, page 51] and 
therefore increase in the value of Kg does not decrease the 
selectivity ratio appreciably ’ [Fig. 18, 19 and 20 and Appendix 
A- 18, A-19 and A- 20], 

The curves for K^/K^ values of zero are close to those 
for K-/K. values of 0,25. It is interesting to note that the 

«*• ' . .TV 

rate equations reduce to simple first and second order when 






Fig. 19 -Selectivity ratio,SR,cis a function of modified 
Thiele parameter,®*^, for selected values of rWK, 



84 


0 » Series R eactions 

The discussion is based on reactions 
A ^ 13 -£• G 
33 '> 3D + 0 

and the selectivity ratio is defined as the ratio of selectivity 
with and without diffusion effects. The selectivity is defined 
as the ratio of the rate of formation of B to that of disappear me j 
of A, that is, r B/" r A * Bor the reaction system considered, the 
selectivity without diffusion limitation is unaffected by any 
of the variable discussed below, since the product concentration 
at the pore month is zero. 

When the first reaction is highly reversible corresponding 
to lower values of equilibrium constant K^, selectivity ratj r 
increased as a result of diffusion limitation. The rate of 


formation of product D depends upon the concentration of inter- 
mediate component B and as first reaction is more reversible, 
the rate of increase of concentration of component B is reduced 
in diffusion limited region. The net effect of decreasing the 
value of at a given value of modified Thiele parameter is 
to reduce the rate of formation of component 3D much more than 
that of component B. Therefore the selectivity with diffusion 
limitations is increased which corresponds to higher values of 


selectivii 


II 1 H 


o [Big. 21 and Appendix 


A- 21] . 
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When the second reaction is faster, than the first, 
corresponding to Ion values of the intermediate product 33 
is converted to product 33 at a 'much rapid rate and therefore 
the selectivity with diffusion limitation is decreased wMch 
corresponds to low values of selectivity ratio [Fig. 22 and 
Appendix A-22] . 

Fig. 23 [Appendix A-23] gives the value of selectivity 
ratio as a function of modified fhiele parameter for selected 
values of K^p o . Hie selectivity ratio is increased with increase 
in the value of i^p Q for any value of 0^. The effect of varying 
K^p q becomes more apparent if (l+K^p Q ) is talcen out of the 
modified Thiele parameter and the values of selectivity ratio 
is presented as a function of ordinary Thiele parameters Whe~ 
presented this way for selected values of VIII, tir 

spread between the curves is increased quite significantly * 


Table V III 

Selectivity Ratio, S R , as a Function of Ordinary 
Thiele Parameter, for Selected Values of K^p Q 










' ' ' . : • . - 

Fig. 23 -Selectivity ratio, Sr, as a function of modified 

Thiele parameter,(p Mj for selected values of KaPo* 


ic 


0-6 

cn 
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Diffusivity of component B is most important in series 
reactions. If diffusivity of B is large, component B will 
diffuse out of the pore at a faster rate and therefore, its 
concentration inside the pore will decrease. This will increase 
the formation of B while reaction of 33 to product D will decrease 
therefore, increase in the value of 33^/33^, for given value of 0 ^ , 
will increase selectivity ratio quite significantly. Conversely 
the decrease in the value of Bg/B^ will decrease Sg significantly 
[Big. 24 and Appendix A-24]. 


When diffusivity of component C is increased, its concen- 
tration inside the pore is decreased resulting in an increase in 
the net rate of first reaction, and consequently increase in the 
concentration of component B. Higher concentration of component. 

B results in increased conversion (irreversibly) to produce D 
and C according to the second reaction, the formation of componei 
G as per second reaction may partly compensate for the decrease 
in the concentration of component C due to the increased diffu- 
sivity. The rates of formation of components B and 33, that is 
r-g and r^ are increased in absolute terms but their ratio 
is decreased with increase in B^/D^ at any given value of 0^. 

The net result is, therefore, to decrease the selectivity ratio 
with increase in the diffusivity of component C [Big. 25 and 
Appendix A-25]* 
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The effect gf change in the value of diffusivity of 
component D irj similar to that for independent and parallel 
reactions, and ia small in comparison to components B and G. 
Increase in the diffusivity of product 3) decreases its partial 
pressure and this decreases the value of denominator of the 
rate expressions and therefore decreases the reaction rate for 
Both the reactions. The effect of diffusion is thus increased, 
which reduces the selectivity ratio as the diffusivity of 
product D is increased [Pig. 26 and Appendix A- 26], 

The effect of varying the adsorption equilibrium constant 
for component B, C or D is similar to that for independent 
reactions. Increase in the value of adsorption equilibrium 
constant increases the selectivity ratio for a fixed value of. 
modified Thiele parameter. The effect for component G is more 
than that for component B or 1, since Uq is multiplied by the 
partial pressure of component C which is much larger than the 
partial pressure for component B or D as product G is formed 
by both the reactions [Pig. 27, 28 and 29 and Appendix A-27, 

A- 28 and A- 29]. 





U) 



Fig. 27-Select 


Thiele para i 



Fig. 29- Selectivity ratio, S r, as a function of modified 

Thiele parameter, 0 m . for selected values of Kq/Ka. 

If ; - I ■■ a : V ; : 
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CH APT3R VI 


C01TCLUSI0HS 


and 


IffiCOIMENDATIONS 


A. Co nclusio ns 


A direct netiod to obtain the missing initial conditions 
of a system of nonlinear differential equations with two point 
boundary conditions using invariant imbedding concept is developed 
which involves less computational error and time. She only 
limitation of t: is new technique is its inability to use very 
small increments, say, less than 0.1. The technique is used to 
obtain the value of overall reaction rates of different products 
in a pore for complex reaction systems. The technique gives 
satisfactory results for modified Thiele parameter values upto 
two. For higher values of 0^, this method does not give accurate 
results as the increment used is not sufficiently small. Quasi- 
linearization technique is used to solve these equations for 
values of 0L upto a value of ten. However this technique has 
also failed in many cases, especially for series reactions, 
for 0^ values larger than three. 


The selectivity of a complex reaction is significantly 
affected due to diffusion limitation ina catalyst pore. The 
magnitude depends upon the type of the reactions, reaction 
mechanism and the physical and chemical constants influencing 

the system. The effect of different variables investigated in 
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this work oil the selectivity ratio is as follows; 

i. The selectivity ratio is found to decrease, 
except for two cases, with increase in modified Thiele 
parameter up to a value of five. Further increase in the 
value of pEjnj does not effect the selectivity ratio. The 
selectivity ratio is found to increase with increase in 
0^ only for independent reactions './her. the value of variable 
D-g/D^ or Kg are less than unity, 

ii. For values of modified Thiele parameter less than 
0.2, diffusion limitations are insignificant and selectivity 
ratio is closo to unity, 

iii. The selectivity ratio is found to increase with 
increase in the value of for parallel as well as indepen- 
dent reactions. Contrary to this, selectivity ratio decreases 
with increase in the value of K-^fox series reactions. However 
the effect of change in IC^ is significant only for 0^ values 
greater than unity. The effect of thousand fold increase 
in the value of K, has maximum effect for the independent 
reactions and minimum for series reactions. In case of 
independent reactions, for value of 100, that is, first 
reaction is also essentially irreversible, the value of 
selectivity ratio remains close to unity even for 0^ values 
as large as ten. 



99 


iv. Mien the second reaction is faster than the first, 
that is, values of less than unity, the selectivity ratio 
for independent reactions is always greater than unity and 
is found to increase with the increase in the value of ^ or 
decrease in the value of K^. For the value of K ^ greater 
than unity, selectivity ratio rapidly decreases with increase 
in the value of ^ or In case of parallel and series 

reactions, lower values of selectivity ratio are obtained 
as the value of decreases from 10. to 0.1 and $2^ value 
increases. 


v. Increase in the value of K.p decreases the effect 

A o 

due to diffusion for all types of reactions studied, but the 
effect in each case is quite small for a given value of 0^. 

For very small value of K^p Q and for the range of K^/K^ values 
investigated, the rate expressions becomes simple first 
, or second order and $2L values are close to ordinary Thiele 
parameter, 

vi. The effect of increase in the value of diffusivity 
ratio Dg/D^ is significant for all the three reaction 
schemes. The selectivity ratio increases with increase in 

: nh h;':n -in- h; : y ■ : n : : :■ . 

the value of Ih/jD. . The maximum effect is observed for 
J5 A 

series reactions while the-, effect for parallel reactions is 
the least for all values of exceeding 0.5. 
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vii. She effect of increasing the diffusivity ratio 
Dq/D^ is to increase the selectivity ratio for independent 
and parallel reactions hut to decrease it for series reactions. 
For values of 0L less than unity, parallel reactions have 
largest effect followed hy independent and then series 
reactions. However for large values of 0^, maximum effect 
is shown hy independent reactions followed hy series and 
then parallel reactions. 

viii. She effect of increasing the diffhsivity ratio D^/D^ 
is comparatively insignificant for all the three reaction 
systems. Increase in the value of D^/D^ decreases the 
selectivity ratio for independent and series reactions while 
its value is increased for parallel reactions. The effect 
increases with increase in the value of 0^ for the first 
two reaction systems, hut in case of parallel reactions, 
maximum effect is observed at value close to unity. 


ix. The effect of increase in the diffusivity ratio 
of the two reactants Bg/D^ in case of independent reaction 
is to decrease the selectivity ratio. The value of selecti 
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x. She eiiect of sixteenfold increase in the value of 
adsorption equilibrium constant ratio ICg/K^ does not have 
any significant effect on selectivity ratio. With increase 
in the value of the selectivity ratio is found to 

increase slightly for independent and series reactions with 
a maximum effect at ^ value close to two* The selectivity 
ratio is found to decrease slightly for parallel reactions, 
with increase in the value of Kg/K^ with a maximum effect o 
at $L close to 0.7. 

xi. The effect of varying the adsorption equilibrium 

is essentially similar to 

that of Kg/K^ hut the magnitude of effect is comparatively 
more for followed by and then Kg/K^. 

xii. The effect of the change in the adsorption equ : ". ’ ■ ... 
oonstant for component S in the oaoe of independent reactions 


constant ratios K^/K^ and Kg/E 


is somewhat peculiar. The selectivity ratio increases with 
increase in the value of Kg/K^ if the value of ^^Less than 
two, but for values of larger than 2.3, Sg increases with 


decrease in the value 


of k b A a . 


However the effect is 


small in magnitude 
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B. Re c ommenda ti oas 

She computed values of selectivity ratio are calculated 
by varying only one variable at a time while keeping the values 
of all other variables at unity# As such these curves may not 
be applicable to any actual system and it is therefore necessary, 
that a method is evolved to estimate the effect of diffusion on 
selectivity ratio when large number of variables have values 
other than unity. Once such a procedure is evolved, then any 
actual system can be handled using only limited data as presented 
in this thesis. A further extension of the work can also include 
the study of the effect on selectivity ratio of different concen- 
trations of reactants and products at the pore mouth and then 
it will be possible to estimate the overall selectivity ratio 
for an integral reactor* 
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APPEIIDXX A-l 

Selectivity Patio, S^, as a function of Modified 
Thiele Parameter, $2^ , for Various Values of 


Reaction A^=±B + 0 
E — + 0 
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AP PENDIX 1-2 

Selectivity Patio, S^ f as a Function of Modified 
Tliiele Parameter, for Various Values of Kg 
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Selectivity Ratio, S^, as a Function of Modified 
Thiele Parameter, 0L, fox Various Values of K,n 
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Selectivity ilatio,S R , as a Function of Modified 
Thiele Parameter, for Various Values of Dg/D^ 
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Selectivity Patio, S^, 
Thiele Parameter, 0^, 


as a function of Modified 


for Various Values 


of B c /B a 


Reaction: A 7=hS + 0 
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1.000 

1.000 

1.000 

.4 

.991 

.995 

.997 

.999 

.999 

.7 

.950 

.969 

.982 

.990 

.995 

1. 

.890 

.925 

.954 

.974 

.986 

2, 

.738 

.802 

.863 

.914 

.951 


.693 

.763 

.834 

.895 

.938 


.6 75 

.751 

.826 

.888 

.934 

10, 

*658 

.748 

.823 

.887 

.933 
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APPENDIX A -6 

Selectivity Patio, S„, as a Punction of Modified 
Thiele Parameter, 0 ^, for Various Values of D^y/jD^ 


Reaction; A-^zz 2 -B + 0 
E — ^D + 0 


\ V d a 

v \ 

.25 

.50 

1. 

2. 

4. 

.1 

! i.ooo 

1.000 

1.000 

1.000 

1.000 

.2 

1.000 

1.000 

1.000 

1.000 

1.000 

.4 

.997 

.997 

.997 

.997 

.997 

.7 

.935 

.983 

.982 

.982 

.981 

1. 

.962 

.957 

.954 

.952 

.951 

2, 

.835 

.872 

.863 

.359 

.856 

3. 

.852 

.841 

.834 

.831 

.829 

4. 

.839 

.831 

.826 

.823 

.821 

10. 

.834 

.828 

.823 

,321 

.820 

- -J 
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APPENDIX A-7 

Selec olvi u/ aatio, S^, as a Etmciion of Modified 

Thiele Parameter } ft for Various Values of D-p/D. 

^ A A 

Reaction; B + 0 

E — =>D + G 


' S '\^e/ D A 

\ x. 

.25 

* 50 

1. 

2» 

4. 

.1 

1.010 

1.004 

1.000 

.998 

.997 

.2 

1.039 

1.013 

1.000 

.994 

.990 

.4 

1.143 

1.045 

.997 

.9 73 

.962 

.7 

1.346 

1.104 

.982 

.922 

.892 

1. 

1.517 

1.148 

.954 

.857 

.809 

2. 

1-709 

1.190 

.863 

.683 

.593 

3, 

1.728 

1.192 

.834 

.613 

.497 

4. 

1.730 

1.193 

.826 

.589 

.454 

10. 

1.735 

1.194 

.823 

.574 

.412 
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APPENHEX A- 8 

Selectivity Ratio, S^, as a Function of Modified 
fhiele Parameter, 0^, for Various Values of Kg/K^ 

Reactions A + 0 

E > D + G 


K. w 1 ^ 

r\ 

.25 

.5 

1, 

2, 

4. 

,1 

1,000 

1.000 

1,000 

1,000 

1.000 

, 2 

1.000 

1.000 

1,000 

1,000 

1,000 

.4 

.997 

.997 

.937 

.997 

,997 

.7 

.981 

.982 

,982 

.983. 

,985 

1. 

.951 

.952 

>954 

.957 

.962 

2.- 

.853 

,860 

.863 

.870 

.881 

3. 

.831 

.832 

.834 

.839 

.847 

4. • 

.823 

.824 

.826 

.829 

,835 

M 

; o 

• 

.821 

.822 

.823 

.826 

.831 
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APPENDIX A-9 

Selectivity- ilatio, S^, as a Function of Modified 
Thiele Pursue ter, 0^, for Various Values of K^/K^ 

Reactions A 3 + 0 
E — 3 + 0 


N. > 

\ 

4 \ 

! 

0.0 

.25 

.50 

1. 

2. 

4. 

.1 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

* 2 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

.4 

.937 

.997 

.997 

.997 

.997 

.998 

.7 

.900 

.980 

.981 

.982 

.984 

.987 

1. 

.945 

.948 

.950 

.954 

.960 

.968 

2. 

.846 

.851 

.855 

.863 

.877 

.898 

3. 

.822 

.825 

.828 

.834 

.845 

.861 

4. 

.016 

.819 

.82! 

.826 

.834 

.847 

10. 

.815 

.817 

.819 

.823 

.830 

.839 
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APPEiTDI X A— 1 0 

Selectivity Eatio, S R , as a Function of Modified 
Thiele Parameter, for Various Values of ¥Ljy/K^ 


Reaction: A p=±.B + C 
E — V D + 0 
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APPEITDIX A- ll 

Selectivity Ratio, S-^, as a Pune ti on of Modified 
Thiele Parameter, for Various Values of K-^/K^ 

Reaction: A ’rnr 1 3+0 

E — » D + 0 


\ 

> 

\ 

l \ 

. - - - - * \ 

.25 

.50 

1. 

2. 

4-* 1 

■ 

1 

. i 

■ ■ \ 
i 


.1 

j 1.000 

1.000 

1.000 

1.000 

1.000 


. 2 

1.000 

i 

1.000 

1.000 

1.000 

1.000 


.4 

.396 

.997 

t 

.997 

.998 

.999 


.7 

i .978 

.979 

.982 

.986 

.992 

1 

M 

! .948 

.950 

.954 

.961 

.973 

2 

m 

.364 

t<~\ 

UO 

CO 

* 

.863 

.865 

.872 

3 

# 

.859 

.837 

.834 

.031 

.829 

4 


.831 

.829 

.826 

.821 

.813 

10 

* 

.829 

.827 

.823 

.818 

.809 
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APPENDIX A-12 

Selectivity Ratio, S^, as a function of Modified 
Thiele Parameter, j2^, for Various Values of 


Reaction: 2A^=^3 + C 

A — >2 + o 


\ K 1 

4 \ 

.1 

.2 

*5 

1. 

2. 

5. 

10. 

100. 

.1. 

.996 

.996 

.996 

.996 

.996 

.996 

.996 

.996 

.2 

.983 

.984 

.984 

.984 

.984 

.984 

.984 

-.984 

.4 

.936 

.940 

.943 

.944 

.944 

.945 

.945 

.945 

.7 

.824 

.847 

.861 

.866 

.869 

.870 

.871 

.871 

1. 

.706 

.748 

.779 

.791 

.798 

.802 

.803 

.804 

2. 

.496 

.543 

.599 

.631 

.651 

.665 

.669 

.674 

3. 

.452 

.488 

.535 

.566 

.592 

.613 

.622 

.629 


.443 

.476 

.516 

.544 

.569 

.593 

.605 

.617 

10 . 

CO 1 

KV 

^ 3 

• ^ 

.470 

.506 

.529 

— 

— 

— 

— 
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APPENDIX A-14 

Selectivity Ratio , S^, as a Function of Modified 
Thiele Parameter, for Various Values of K,.;p 

Reactions 2A B + C 
A — >D +• G 
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APPENDI X A-15 

Selectivity Ratio, S R , as a junction of Modified 
Thiele Paraaeter, 0^, for Various Values of Kg 


Reaction? . 2A^=i 33 + 0 
A — > 33 + C 



.1 

.956 

.978 

.991 

.996 

.998 

.999 

.2 

.357 

.920 

.967 

.984 

.993 

.997 

.4 

■ i 

.670 

.778 

.892 

.944 

.972 

.989 

.7 

■ ' 

r nrr 

.5oo 

.622 

.770 

.866 

.929 

.972 

1. 

.482 

.539 

.676 

.791 

.832 

.950 

2. 

.456 

.472 

.534 

.631 

.752 

.879 

3* 

.449 

.461 

.502 

.566 

.671 

.819 

4. 

.446 

.457 

.492 

.544 

.629 

.772 

10. 

.438 

.451 

.433 

.529 

.594 

.697 



10 . 


.000 

.997 

.995 

.936 

.975 

.936 

.898 

.866 

.777 
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..APPENDIX A-15 

Selectivity Ratio, S^, as a Function of Modified 
Thiele Parameter, 0^, for Various Values of Dg/D^ 



Re ac ti on 

: 2 A ^5 

A - 

f=^b + 

— >• D + 

C 

0 








% \7 

i \ 

.25 

.50 

1 . 

2. 

4. 

.1 

.995 

.996 

.996 

.996 

.996 

.2 

.980 

.983 

.984 

.985 

.985 

.4 

#9 31 

.939 

.944 

.946 

.947 

.7 

.840 

.857 

.866 

.871 

.874 

1 . 

.752 ' 

.776 

.791 

.800 

.805 

2. 

.559 

.601 

.631 

.652 

.664 

3. 

.495 

.533 

.566 

.597 

.614 

4. 

.481 

.514 

.544 

.570 

- 

lo.: 

.473 

.504 

.529 


— 
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APPENDIX A-16 


Selectivity Ratio, S R , as a Function of Modified 
Thiele Parameter, for Various Values of Dq/D^ 


Reaction: 2A + C 

A — *■ D + G 
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APPENDIX A-1,7 

Selectivity Ratio, S^, as a Function of Modified 
Thiele Parameter, 0^, for Various Values of Dp/D^ 


Reaction: 2A^s=*B + C 

A — f D + G 


V d a 

NSVS V 

i \ 

.25 

.50 

r. 

2. 

4. 

.1 

.994 

.995 

.996 

.996 

.997 

.2 

.976 

.981 

.984 

.986 

.986 

• 4 

.920 

.935 

.944 

.948 

.951 

.7 

.830 

.853 

.866 

.873 

.877 

1 . 

.755 

, 778 

.791 

.798 

.802 

2. 

.613 

.626 

.631 

.631 

.633 

3. 

.555 

.567 

.566 

.567 

.564 


.531 

.542 

.544 

.545 

.543 

10. 

.503 

.522 

.529 


— 
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APPENDIX A-18 

Selectivity Ratio, S R , as a Function of Modified 
Thiele Parameter, 0 ^, for Various Values of Kg/K^ 
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AP PENDIX A-,19 

Selectivity Patio, S^, as a Function. of Modified 
fiiiele Pax one ter , 0 ^, fox Various Values of K^/K^ 


lie actions 2A + C 

A - — *-D + C 


V' VtT 







\ ~ ^ 

% \ 
l \ 

0.0 

.25 

.50 

1. 

2. 

.4. 

.1 

.937 

.997 

.997 

.936 

.995 

.992 

.2 

.939 

.937 

.986 

.934 

.980 

.972 

.4 

.957 

.954 

.950 

.944 

.932 

.912 

.7 

oo n 
♦ uu ( 

.302 

.376 

.365 

.849 

.821 

1. 

.811 

.806 

.801 

.791 

.775 

.747 

2. 

.629 

.631 

.632 

.631 

.623 

.608 

3. 

.562 

.564 

.565 

.565 

.565 

.554 

4, 

- 


.546 

.544 

.542 

.526 

10. 

i. 

- 

- 

- 

.529 

.519 

,498 
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APPBIIDIX A- 20 

Selectivity Hatio, as a Function of Modified 
Thiele Paraneier , f° x Various Values of 


lie a c lion; 2k~z=±:'£> + 0 

A *1 + C 


&U 'V 

l \ 

0.0 

.25 

.50 

1. 

2. 

4. 

.1 

.997 

.997 

.996 

.996 

.995 

.994 

.2 

.987 

.936 

.936 

.984 

.981 

.976 

.4 

.953 

.951 

.948 

.944 

.935 

.920 

.7 

o 0*1 
*OO.L 

.377 

.873 

.866 

.853 

.830 

i. ; 

.306 

.302 

.793 

.791 

.778 

.755 j 

2. 

.629 

.630 

.631 

.631 

.630 

.613 

3. 

.552 

.563 

.564 

• 566 

.567 

.557 

1 

4. 

- 

.544 

.544 

.544 

.542 

.533 

10. 

- 

- 


.529 

.521 

1 
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APPENDIX A- 21 

Selectivity Ratio, as a Function of Modified 

Ehiele Parameter, 0^, for Variotis Yalues of K-^ 


Reaction; A- ^= £ B + C 
B — + D + 0 


0JI 

L 

Kl 

' . \ 

.2 

.5 

1. 

2. 

5. 

10. 

100. 


.1 

.997 

.997 

.997 

.997 

.997 

.997 

.997 

.997 


.2 

.987 

.987 

.987 

.987 

.987 

.987 

.987 

.987 


.4 

.952 

f . 

.952 

.952 

.952 

.952 

.952 

.952 

.952 


.7 

.882 

*880 

.879 

.879 

*878 

.878 

.878 

.878 


1. 

.816 

.810 

.806 

.804 

.803 

.802 

.802 

.802 


2. 

.692 

.675 

.655 

.646 

.640 

.636 

.635 

.635 


3. 

.64 

.624 

.600 

.588 

.580 

.573 

.570 

.570 


4. 

V 

.60 

- 

.563 

- 

— 

— 

- 
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APPEH DIX A- 2 2 

Selectivity' Ratio, S^, as a Function of Modified 
Ihiele Parana ter, fox Various Values of 


Reaction: A ^=1 3 + 0 
B — * D + C 



A .687 .806 .910 .952 .975 .990 .995 

.7 .484 .624 .790 .879 .934 .972 .986 

1. .338 .511 .638 .304 .087 .950 .974 

2. - .29 .381 .526 .645 .760 .875 .930 

3. I - - .583 .696 .823 .893 

4. - ! - - - .66? .791 .867 

I i 
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APPENDIX A-23 

Selectivity Ratio , S^., as a Function of Modified 
Thiele Parameter, 0^ f for Various Values of 


Pteaction; A -^±33 + 0 
B ♦D + C 
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APPENDIX A- 24 


Selectivity Ratio, S E< as a Pune ti on of Modified 
Thiele Parameter, 0^., for Various Values of Dg/D^ 


Reaction; + C 

B *3) + C 
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. .-A PPENDIX A- 2 5 


Selectivity Ratio, S R , as a Function of Modified 
Thiele Parameter, 0^, for Various Values of 


Reaction: A B + 0 
B — *D + C 
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A PPEK DIX A- 26 

Selectivity Patio, S^, as a Function of Modified 


•/jOliiele 

Parameter , 

0 M > for 

Various 

Values 

of vy 


Reactions A^3 + C 





B - 

D + C 



V“T" 

| 

i 

| '■ * . 





V* 

V,,\ 

| .25 

j 

.50 

1 . 

n 

<£. « 

4 . 

4 X % j 






.1 

! .997 

.997 

.997 

.997 

.997 

.2 

.987 

.937 

.987 

.937 

.987 

.4 

i .952 

.952 

.952 

.952 

.952 

.7 

.382 

.880 . 

.379 

.378 

.878 


; .314 

.808 

.804 

.802 

.801 

■ i 

. i 

2. 

j . .677 

.658 

. 646 

.639 

.635 

3. 

| .523 

i 

.602 

.538 

.530 

.576 1 

■ ! 

4. 

f ■ : - • ; . ; 

! .593 

- 

~ - 

- 

i 


i 
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AP PEND IX A-,2 9, 

Selectivity Ratio, S^, as a Function of Modified 
Thiele Parameter, fox Various Values of I^/K^ 


Reactions A^nS + C 
B — * D + C 
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appendix b-i 


COMPUTER PROGRAMME FOR INVARIANT IMBEDDING TECHNIQUE 


T JOB CGF04-1,TIME008,NAME A B LAL PAR REAC Kl K2 KAPO 

$ IBJOB 
$ IBFTC 

DIMENSION U(4,4) ,V(4,4) ,W(4,4) ,x( 3 ) »AA(2 ) ,BB(2 ) ,CC (2 > 

READ 61 » ( (U(K»L ) ,L=1 »3) ,K=1*2 ) 

61 FORMAT ( 3 F 1 5 » 8 } 

C0 = 0. 

D0 = 0. 

COMMON A*B1»B2,B3 
MM = 0 

2 READ 51, A,B1*B2,B3 
51 FORMAT ( 4-F7.2 ) 

M=0 

DO 99 J= 8 » 1 0 

N=0 

N1 = 0 

N2=0 

N3 = 0 

N4 = o 

N5 =0 

R = J 

Q=l. /R 
P=Q*Q 

c=co 

D = DO 

5 N=N+1 

Cl=C-P*F (C.DI 
D1=D-P*G(C»D) 

C2 = C1-P*F<C1*D1 ) 

D2=Dl-P*G(Cl,Dl ) 

Cl=C2-p*F(C2,D2 ) 

D1=D2-P*G(C2,D2 ) 

GO TO (7» 28, 7, 29, 7, 1*7. 12, 7, 3, 7, 3, 7, 1,7, 14, 7, 6, 7, 6, 7, 1,7, 6, 7, 3, 7, 3, 

1 7, 1,7, 21, 7, 8, 7, 8, 7, 1,7, 8, 7, 3, 7, 3, 7, 1,7, 8, 7, 6, 7, 6, 7, 1,7, 6, 7, 3, 7, 3, 

2 7, 1,7, 19, 7, 4, 7, 4, 7, 1,7, 4, 7, 3, 7, 3, 7, 1,7, 4, 7, 6, 7, 6, 7, 1,7, 6 ,7, 3, 7, 3, 

3 7, 1,7, 4, 7, 8, 7, 8, 7, 1,7, 8, 7, 3, 7, 3, 7, 1,7, 8, 7, 6, 7, 6, 7, 1,7, 6, 7, 3, 7, 3, 

4 7, 1,7, 25), N 

7 C3=C-P*(F(C»D)+F<C1,D1)+F(C2,Q2) ) 

D3=.D-p*(G(C,D)+G(Cl >Dl ) -t-G ( C2 »D2 ) ) 

C = C3 
D = D3 
GO TO 5 
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28 C=CO 
D = DO 

11 C4'=C-P*( F(C,D)+F ( Cl »01 )+F (C2»D2 )+F( C3»D3 ) ) 

D4 = D-p*(G'(C,D)+G(Cl,Dl ) +G (C2 ,02 ) +G( C3»D3 ) ) 

C = C4 

D = D4 
GO TO 5 

29 C=C0 
D = D0 

15 C5 = C-p*{F(C,D)+F(Cl,Dl )+F(C202 >+F{C3,D3>+F(C4,D4) ) 
D5=0-p*(G(C,D) 6{C1,D1 ) +G { C2 , 02 ) +G( C3 ,03 ) +G (C4 , D4 ) ) 

1 C = C5 
D = D5 
N1=N1+1 

60 TO(5»ll»5»ll»5*ll»5»ll*5»31*5»ll»5»ll»5*ll» 

1 5»11»5»11»5»11»5»11*5»31»5»11»5»11»5*11)»N1 

12 C=CO 
D=DO 

20 C6=C-p* { F { C,D ) F(Cl,Dl)+F(C2»D2)+F{C3,D3)+F(C4»D4)+F(C5,D5n 
D6=D-p*(G(C,D) +G ( Cl »D1 ) +G( C2 , D2 ) +G( C3,D3 ) +G ( C4,D4 )+G( C5 ,D5 > ) 

3 C=C6 
D=D6 
N2=N2+1 

60 TO (5, 11, 15, 5, 11, 15, 5, 11, If, 5*11, 15, 

1 5 ,11 ,15,5, 11, 15,5, 11 ,15, 5, 11 » 15 > »N2 

14 OCO 
D = D0 

18 C7 = C-p* ( F { C ,D ) +F ( Cl , D1 ) +F ( C2 » D2 > + F ( C3 ,03 ) +F ( C4 , D4 ) +F ( C5 ,D5 ) + 
1 F(C6,D6 ) ) 

D7=D-P*(G(C,D)+G(C1,D1 )+G(C2,D2»+G(C3,D3 > +G ( C4 » D4 ) +G ( C5 , D5 > + 
1 G(C6,D6 ) ) 

6 C=C7 
D=D7 
N3=N3+1 

GO TO (5, 11, 15, 20, 5, 11, 15, 20, 5, 11, 15, 20, 5, 11, 15, 20 ) ,N3 

21 C=C0 

D = D0 ■ ■ \ 

22 C8=C-P*(F(C,D)+F (C1»D1)+F(C2,02 )+F(C3,D3 )+F{C4,D4)+F(C5 ,D5)+ 
1 F( C6,D6 >+F(C7,D7 ) ) 

D8 = Dr-p* < G ( C ,D ) +G CCZ , D1 } +G ( C2 , D2 > +G { C3 , D3 > +G { C4 , D4 » +G ( C5 , D5 ) + 
1 G(C6,D6)+G(C7,D7 ) ) 

IF ( J-8 )8 ,40,8 
40 V( 1,4)=(C8-C0)/Q 
W( 1,4)=(D8-D0)/Q 
GO TO 99 
8 C=C8 
D = D8 
N4=N4+1 

GO TO (5, 11, 15, 20, 18, 5, 11, 15, 20, 18 ) ,N4 
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.9 C = C 0 

C9=C-p*(F(C,DI+F(Cl,[)l)+F(C2,D2l+F(C3,D3)+F(C4,D4l+F(C5,D5l+ 

1 D F 9 = D-;“s!c;dUgIcUD1);GK2,02.FG(C3,D3, + 3<C4,D4, +S( C 5 ,D3,f 

1 6 ( C6 »D6 ) +G ( C7 » D7 ) +G ( C8 » D8 ) > 

I F ( J-9 ) 4 , 41 ,4 
+1 V(2,4)=(C9-C0)/Q 

W ( 2 » 4 ) = ( D9-D0 ) /Q 
GO TO 99 
+ C-C.9 

D = D9 
N5=N5+1 

25 ?; 0 IS^“fi?o'D5^«:0?)*F,C2.D2. + F(C3.D3.*F,C4.D4) + F,CS.D 5 » + 

i d ^^^!S:^!:^^^::^^:o^g,c3,o3 1+ g,c4,o4, + g,c3, D 3.f 
1 G ( C6 »D6 ) +G ( 07 »D7 ) +G ( C8 »D8 >+C ( C9 » D9 ) ) 

V(3,4)=(C10-C0)/Q 
W(3,4)=(D10-D0)/Q 
99 CONTINUE 
68 DO 62 K= 1 » 3 
DO 62 L= 1 » 3 

62 V(K,U=U(K,L) 

DO 63 1=1,2 
11=1+1 

DO 63 K= I I ,3 

63 V(k!l)=V(k!l)-V(I,L!*V(K,I)/V(M) 

DO 64 LQ=3,4 

LL=6— LQ 

L_ J = L- L — 1 ^ 

6 4 V^K^4 ) =\M >-V (LL,4)*V(K»LL) /V ( LL » LL ) 

DO 65 K = 1 » 3 

65 X(K)=V(K,4)/V(K-*KI 
DO 67 K= 1 , 3 

67 V ( K»4) =W ( K,4) 

M = M+ 1 

A A ( M ) = X ( 1 ) 

' - ■ , fe'BT 

' ' : , ,’CC t Hi *XJTS> * 

IF(M-2)68,98,98 )/((? * AA ( 1 > +AA f 2 ) ) *B2*B3 J 

96 suMi«;^i;«, l ..c c « i ..B B «2..cc,2 1 

69 FORMAT ( 4F7.2»Fl4,4,Fl6«6 ,Fll * 3 » 

MM=MM+1 

I F ( MM-200 ) 2,97 ,97 

97 STOP 
END 
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$ I BFTC FUN1 

FUNCTION F(C*D) 
COMMON A»Bl*B2,B3 
F=-A#A*B3/(1.+B3>*( < 


1.-2.*CtR>**2-C*<C+D>/B1)/U.+B3*D/(1. + B3>>**2 


return 

END 

SIBFTC FUN2 

FUNCTION G ( C » D ) 


COMMON A*B1jB2,B3 . R - 

G=-A*A/B2*(1.-2.*C-D)/(1.+B3*D/(1.+B3 

RETURN 

END 



sentry 

1 . 

1 . 

1 . 


.125 

.11111111 

.1 


.015625 

.01234568 

.01 



139 


APPENDIX B-2 

COMPUTER PROGRAMME EOR QUAS I L IN EAR I ZAT ION TECHNIQUE 


$JOB CGF0A1,TIME008,NAME A B LAL IND REAC Kl K2 KAPO 

$ I B JOB 

$IBFTC dimension X , 1000 ). Y.1000).H(1000..E(I0001.P(1000.. 0(10001.6.1000. 
1 PP ( 1000 ), QQ( 1000), 6G( loop) *R( 1000) ,S( 1000 >,V( 1000 ),WC1000 


NN = 0 

B=. 1E+0T 

11 READ 1,A,B1,B2,B3 
1 FORMAT ( 4F7 . 2 ) 

X0 = 0. 



Y0 = 0. 

B5=B3/(1.+2.*B3) 

DO 2 1=1,1000 
X(I)=O.OOl 
Y(I) =0.001 


DO 3 1=1,999 
R( I )=X( I )-< «5 


— X ( I ) )*(!•- 


X ( I ) — Y { I ) ) / B 1 



S( I )=1.-B5*(X( I )+YM ) ) 

H ( I ) =-( 2 .*B+A*A* ( R ( I ) *B5/S( I )** 2 + ( 1 
F( i )=_(2.*B+A*A/B2*(B5*Y ( I ) /G ( I ^**2+1 
V( I )=A*A*(R( I )/S( I )~X( I )*(., ' 

/B1 ) / S ( I )■) ) 

W ( I )=A*A/B2*(Y( I ) / S C I ) 

W(1)=W(1)-.5*B 
V(l)=V(i)-.5*B 
H(999)=H(999)+B 

E ( 9 99 ) = E ( 9 99 } +B 
P( 1 )=H(1 ) 

Q(1)=B/P(1) 

G( 1 )=V( 1 ) /P ( 1 ) 

PP(1)=E(1) 

QQ ( 1)=B/PP( 1) 

GG( l)=Wt D/PPt 1 ) 

^ ~ a * — O OOQ 


+( 1*5-2#*X( I )-YC I } 

- .7 s ( i > ) ) 

(R(I)*B5/S( I)**2+(l. + U.5-2.* 
Y ( I )*<85*Yf I )/S( I ) **2+l*/S( 1 ) ) ) 


! — 1> ) /P( I ) 

- 

G(I-l) )/PP<D 


) 
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DO 5 1=1,998 
M=999— I 

X(M)=G(M)-Q(M)*X (M+l ) 

5 Y(M)=GG(M)-QQ(M)*Y(M+1) 
Xl'=(X(l)-.5)*1000. 
Y1=(Y(1)-.5)*1000. 

IF { ABS(XO-Xl )-.001 ) 15,15,16 

15 IF ( ABS ( Y 0-Y 1 ) - • 001 ) 17,17,16 

16 XO=X 1 
YO = Y 1 
GO TO 7 

17 AB=X1/(Y1*B2) 
PRlNT6,A,Bl,B2,B3, AB,Xl,Yl 

6 FORMAT (4F7.2,3F15.5 ) 

8 NN=NN+1 

I F ( NN-99 ) 11,88,88 
83 STOP 
END 

SENTRY 
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APPE17DIX G 


ILLUSTRATIV E KJ&PLE: DEHYDRATI OIT OE ETHANOL 


Butt et.al [39] considered the dehydration of ethyl 
alcohol on alumina catalyst and obtained the following rate 
expressions at 274°C in gram moles per min, per gram of catalyst. 

0.00627 x 0.00135 p! 

B 4.x (1+0.00135 P A + 0.001 p B + 0.00151 P 0 ) 2 

( la) 

0.000099 x 0.00135 p A . 

r D. = XT+oZooi 33 ~v A +~o 7 obTp B To'.' 00151' pj) 

where the reactions can be expressed as 

2A ?=** 13 + C 
A * D + 0 


Component A stands for ethyl alcohol, B for diethyl ether, G for 
water and D for ethylene. The partial pressures in the above rate 
expressions are in mm, of Hg. The density of alumina catalyst 
pallet is 0.80 gram/cc and porosity is 0.667. Changing the 
partial pressure in atmosphere and rate of reactions in gram 
mole per sec. per cc of pore volume, the following rate expressions 


are obtained: 


an 


0,000033 P 


r B " 


A 


(l + 1.026 p A + 0 .7 6 Pj + 1*1576 p Q ) 


(2a) 
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0,00000203 


^ (1 + 1.026 p A + 0.-76 p B + 1.1576 p 0 ) 2 


(2b) 


Here 


= 1.026 
= 0.76 
= 1.157 


= 0.101 
= 0.0736 
= 0.216 
= 0.139 


Therefore K a p q = 1.026 

1c l E a = ,000033/1.026 = .0000322 


K i = 00 


K 2 . = ^A 2 = ^ggl§3 - 15.86 

1+K. p = 2.026 

A o 

In this problem, since the values of all the variables 
are different than unity , the figures given earlier cannot be 
used directly and the problem has to be solved by Computer. 
Equations 25a and 25b of Chapter II can be written as 


*The values for effective diffusion coefficients are taken as 
one-half of the bulk diffusion coe^-xicients. 
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dr) 


?B = 
2 


(a 

K A p o ‘ 1 h 


k^Ri 1 


D^ 


( 1 + Vb + 6 1 9 D^' 


(3a) 


£^__ T a ¥aL 

dr] 2 o ) 


2D* Dp, 

^ a l “ 9 B “ fj C?I) 
i + Vb + \ C ?P 


(51) 


where a-p = 1.0 

' k aP„ 
A^o 


_ . , h _ 3. + 5s ^ m 

T i = < ir ; \ k a V 


. ' k a p o / V h fc h ) 

h - iscr; s *i - d a + k a d c 


<c^-| — 1 *0 


-0.331 


(4a) 

(4b) 

(4c) 

(4d) 


Substituting -these values .4h- the above equations, the following 
are obtained. 


d 2 cjXg ■* 4.91 x L 2 (l — 1.46 9 B ” 1*375 ?p) 


diV 


(1 - 0.171 q>B - 0.331 <fy) : 


(5a) 


;£jV 

dt] 2 


..323 x h 2 ( 1 - 1*46 ^ - 1 / 375 _Jsl 

~(l - 0.171 9 b - 0.331 <P D ) 


0.076 


: dc? s 


S B ~ 




T)=0 


2 x 0.076 


dr) 


dc ?D 

r)=0 + 0,159 dr) *" 


T )=0 


5b) 


■ ci-h'; 
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S, 


K 2 K A p o 


1 + K A p n (l + 2 K 2 ) 


= 0,471 


and 


A^o 

S E = s i/ 0 * 471 


Ihe values of selectivity matio is obtained by solving equation 
5a and 5b with the following boundary conditions: 


<P B (°) = 9j)(o) = 0 
T)=l 


dr) 


% 

~ drj 


r)=l 


0 


and for different lengths of the pore, L, She results are 
summarized in the table below: 


is 


1" 

0m. 

i " ' , • 



0,1 

• 226 

.998 

0,2 

.452 

. *995 I 

' ' ^ - ■■■ : , j 

0.5 

; 

1.35 

.980 : 

j 

2.66 

.959 ; 
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